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Abstract. One of the many problems to which J.S. Dowker devoted his attention 
is the effect of a conical singularity in the base manifold on the behavior of the 
quantum fields. In particular, he studied the small-i asymptotic expansion of the 
heat-kernel trace on a cone and its effects on physical quantities, as the Casimir 
energy. In this article we review some peculiar results found in the last decade, 
■ regarding the appearance of non-standard powers of t, and even negative integer 

powers of logi, in this asymptotic expansion for the selfadjoint extensions of some 
symmetric operators with singular coefficients. Similarly, we show that the £- 
function associated to these selfadjoint extensions presents an unusual analytic 
structure. 

1. Introduction 

In Quantum Field Theory, the transition amplitude for particles interacting with 
a background field or subject to boundary conditions can be described in terms 
of the effective action of the model under study. In the one-loop approximation, 
this effective action can be expressed in terms of the functional determinant of the 
differential operator appearing in the quadratic in the quantum field term of the 
action. Since the spectra of these operators are unbounded, these determinants 
must be defined in the framework of an appropriate regularization. 

In 1976, J.S. Dowker and R. Critchley p] presented a powerful and elegant regu- 
larization scheme for the definition of these functional determinants, based on which 
is generically known as the (^-function associated with the differential operator, built 
from the spectrum of the quantum fluctuations. 

Since then, this formalism has been successfully applied to the determination of 
one-loop effective actions, vacuum energies, anomalies and other physical quantities 
of interest in Quantum Field Theory. It is presently a fundamental tool in the study 
of quantum effects in systems under the influence of external conditions [21 El IH El 

To be more specific, the first quantum corrections to the effective action are, 
in general, given by (a power which depends on the nature of the fields of) the 
functional determinant of an elliptic differential operator on an appropriate Hilbert 
space -the quantum fluctuations operator- that is obtained as the second functional 
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derivative of the classical action at vanishing quantum fields. J.S. Dowker and R. 
Critchley gave a definition of this functional determinant in terms of the derivative 
of the associated ^-function, thus relating the one-loop effective action with the 
spectrum of the quantum fluctuations operator and, consequently, with the external 
conditions applied on them. 

This formalism is also related to the so-called proper time regularization through 
the relation between the (^-function and the trace of the corresponding heat-kernel 
[S] (See Appendix |A}. Let H represent an elliptic boundary problem on an m-dimen- 
sional compact manifold M with smooth boundary dM, with a discrete spectrum 
{A n }neN- Then, under appropriate conditions [9], the ^-function is defined as the 
trace of the complex power H~ s , 

(1.1) ( H (s):=TrH- s :=J2 X n, 

rieN 

series that converges to an analytic function in an open half-plane with (s) large 
enough and admits a meromorphic extension to the whole s-plane. On the other 
hand, for positive definite H, the trace of the operator e~ tH is defined as [TO] 

e -tx n 

neN 

for t > 0. It has been established [PJ that for a differential operator H of order 
d with smooth coefficients, the trace of the heat-kernel, Tre~'^, admits a small-t 
asymptotic expansion given by 

oo 

(1.3) Tr e~ tH ~ ^ a n ■ t { - n - rn)/d , 

n=0 

where the Seeley-De Witt (SDW) coefficients a n are integrals on the manifold M 
and its boundary dM of local invariants which depend on the coefficients in H, the 
metric on M and the boundary conditions imposed at dM [10]. Let us remark that 
for a differential operator with smooth coefficients on a base manifold M without 
boundary the SDW coefficients a n vanish for odd n. 

Since the function (h{s) in Eq- (H-ll) is related to the heat-trace Tre~ tH in Eq. 
(11.21) by a Mellin transformation [8], the expansion in Eq. (11.31) implies that (h{s) 
has isolated simple poles at 

in — ti 

(1.4) s n = — , with n = 0,1,2,... 

a 

with residues related to the SDW coefficients by (See Appendix [A} 

(1.5) Res{T(s)( H (s)}\ s=Sn = a n . 

In particular, Res [Cff(s)]| s=Sn = when s n = 0, —1, —2, . . . and Cf/( s ) is analytic in 
a neighborhood of the origin. 



SPECTRAL FUNCTIONS OF NON ESSENTIALLY SELFADJOINT OPERATORS 



3 



The functional determinant of H introduced in [TJ can be defined by 
(1.6) logDetff :=-^GrO0 

where | s _>.o stands for the limit of the analytic continuation to a neighborhood of 
s = 0. 

When employed for the description of the one-loop contribution to the effective 
action, this scheme of regularization leads to local counterterms expressed in terms 
of the SDW coefficients (or, equivalently, in terms of the ^-function). This justifies 
the intense research devoted to the ^-function and heat-kernel methods during the 
last decades in relation with its applications to Quantum Field Theory. 

In this context, one of the problems to which J.S. Dowker devoted his attention 
is the influence that a conical singularity in the manifold has on the quantum field 
behavior [TTJ [121 1131 EEH EES EH EE3 HE]- Particles in a conic singularity, situation 
which can be related to the so-called Calogero models [TH], appear in the study of 
quantum fields in a black- hole background [22], in the presence of cosmic strings (23] 
and in condensed matter, for example. The heat-kernel trace on a conic manifold 
has been studied in detail [20, 21J and shown to admit, as in the regular case, the 
asymptotic expansion described by Eq. (jl.3p plus a possible logarithm. However, 
contrary to the case of a regular background, the SDW coefficients a n do not vanish 
in general for odd n, even in the presence of no (other) boundary. 

In any case, the heat-trace on a cone satisfies the asymptotic expansion in Eq. 
(11.31) only if one assumes that the fields are regular at the singular point (see other 
studies and applications of this topic in (2U 1251 I2H [23 I2S1 I2S])- In fact, if one 
considers other (square integrable) behaviors of the fields at the singularity one 
finds that these models present different spectral properties for which the expansion 
given in Eq. (11.31) no longer holds. 

In this article we will review some results regarding the peculiar properties of the 
spectral functions associated to some differential operators with regular singularities 
in its coefficients. In particular we will focuss on some progress made in the last 
years in the understanding of the heat-kernel trace and ^-function properties of some 
(symmetric but not essentially self adjoint) differential operators which are locally 
homogeneous near the singularity, i.e. with the singular coefficient in the potential 
with the same scaling dimension as the highest derivative in the kinetic term. As we 
will see, an essential aspect of these models is the possibility of imposing (for certain 
range of the parameters in the potential) boundary conditions at the singularity that 
break this scale homogeneity. 

The Laplacian on a manifold with a conical singularity is an example of these 
kind of differential operators. The asymptotic expansion of the heat-kernel trace 
corresponding to this operator has been considered, probably for the first time, 
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by A. Sommerfeld [30] and H.S. Carlslaw |3~I] . It was only in 1980 that it was 
pointed out by C.J. Callias and C.H. Taubes [32J that, for these kind of differential 
operators, the heat-kernel trace small-t asymptotic expansion in terms of powers of 
the form t( n - m )/ d could be ill-defined and conjectured that more general powers of 
t, as well as logt terms, could appear. Indeed, it was proved in the following years 
|33j, for some second order elliptic (essentially) selfadjoint differential operators H, 
the presence of £( n - m )/ 2 logt terms in the small-t asymptotic expansion of the heat- 
kernel at the diagonaQ, e~ tH (x, x), with some distributional coefficients with support 
concentrated at the singularity. 

More recently, E. Mooers [31] studied the selfadjoint extensions of the Laplacian 
acting on differential forms on a manifold with a conical singularity and found that 
the asymptotic expansion of the heat-trace contains powers of t whose exponents 
depend on the deficiency angle of the singularity. 

In similar settings, it was shown in a series of articles [3SJ ESI S3 ESI EH] that, 
for some non essentially selfadjoint locally homogeneous differential operators, the 
small-t asymptotic expansion of the heat-kernel trace presents non-standard powers 
of t, i.e. powers with exponents which are not determined by the order of the 
differential operator and the dimension of the base manifold only, as for the regular 
case (see Eq. (11.31) ). but also depend on the coefficient of the singular term in the 
potential. Consequently, in these models, the presence of a regular singularity in 
the potential term of the differential operator leads to non-standard poles in the 
associated ^-function, which lie at positions that depend on parameters other than 
the order of the operator and the dimension of the manifold. 

Later, K. Kirsten et al. [HH SD 1121 S3] considered a limit case of a symmetric 
second order locally homogeneous differential operator, finding that integer powers 
of logt terms appear in the small-t asymptotic expansion of the heat-kernel trace 
and a logarithmic cut in the corresponding ^-function, so correcting an error in 
Appendix A in [37] . 

It is our aim to review these results concerning the non-standard behavior of 
the small-t asymptotic expansion of the heat-kernel trace and the singularity struc- 
ture of the ^-function corresponding to this kind of locally homogeneous differential 
operators. We will consider the case of some symmetric operators with a regular 
singularity in the potential term with the same scaling dimension as the deriva- 
tive term, which makes them to admit a continuous family of selfadjoint extensions 
(SAE) [44] . As we will see, the (local) scaling homogeneity is in general broken 
in the domains of definition of these SAE, and this fact has consequences on the 
behavior of the associated spectral functions which will be studied in the following. 



By "heat- kernel" we mean the kernel e tH (x,y) of the integral operator e 
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In Section [2] we will consider a Dirac operator D (see Eq. (12.11) ) defined on a space 
of two- component functions &(x) with x taking values in the compact segment [0, 1]. 
We will introduce in the operator D a singular term proportional to 1/x, which has 
the same scaling dimension as the kinetic term d x . As we will see, this simple model 
shows the above mentioned characteristic. 

This first order differential operator is not positive definite and one can not define 
the associated heat-kernel. Instead, we will consider the (^-function as defined in 
(11.11) and show that it presents a non-standard pole structure due to the presence 
of the singular term ~ 1/x in £). In so doing, we will first determine the large- 
|A| asymptotic expansion of the resolvent (D — A) -1 , since the powers of A in this 
expansion determine the position of the poles of the ^-function (see [5], for example). 

We will proceed as follows: Firstly, we will construct the resolvents for two partic- 
ular selfadjoint extensions for which the boundary conditions at the singular point 
x = are invariant under the scaling i ex. The resolvent expansion for these 
particular SAE displays the standard powers of A, leading to the standard poles for 
the ^-function. Secondly, we will show that the resolvent for a general SAE is a 
convex linear combination of these special resolvents with coefficients which depend 
on A. This additional dependence on A leads to the non-standard powers in the 
resolvent asymptotic expansion of a general SAE and, hence, to non-standard poles 
for the associated ^-function. The selfadjoint extensions of D are not, in general, lo- 
cally scale invariant at the singularity in the sense that the conditions the functions 
in its domain satisfy near the origin are not invariant under the scaling x — > ex. 
As c — > they tend to the conditions satisfied by the functions in the domain of 
one of the locally scale invariant SAE, and as c — > oo they tend to the other. The 
dependence of the residues at the anomalous poles on the SAE will also be explained 
by a scaling argument. 

This model describes the central idea of our work. We consider differential op- 
erators H which contain singular terms of the same scaling dimension as the high- 
est derivative term; these operators -which we call locally homogeneous near the 
singularity- admit selfadjoint extensions whose domains are characterized by bound- 
ary conditions at the singularity that break, in general, this local homogeneity. 
This introduces in the definition of the SAE dimensionful parameters which are not 
present in its expression as a differential operator. This determines that the large- 
|A| asymptotic expansion of the resolvent-trace of a general SAE, Tr (H — A) -1 , 
presents non-standard powers of A, with dependence on the coupling in the singular 
term. This, in turn, reflects in the presence of non-standard poles in the complex 
s-plane of the (^-function, TyH~ s , and (for H positive definite) in the presence of 
non-standard powers of t in the small-i asymptotic expansion of the heat-kernel 
trace Tre~ tH . Here, by "non-standard" we mean "not determined by the dimension 
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of the base manifold and the order of the differential operator only" , as happens for 
the regular case. 

Indeed, in Section [3] we consider a second order differential operator on the seg- 
ment [0, 1] with a singular potential term of the form to g(g — l)/x 2 . We show that, 
for \g\ < 1/2, one obtains similar results for the associated ^-function as those de- 
scribed for the first order operator case we consider in Section [2j One finds that there 
are two locally scale invariant SAE for which the (^-function presents the usual poles. 
The resolvent of the general SAE is obtained as a convex linear combination of the 
resolvents of these particular SAE, with a coefficient dependent on the A-parameter. 
This implies the presence of ^-dependent poles in the ^-function with residues which 
also depend on the SAE. On the other hand, the small-t asymptotic expansion of 
the heat-kernel trace of the general SAE presents non-standard (^-dependent powers 
of t with SDW coefficients with dependence on the SAE. 

This second order differential operator taken for g = 1/2 also admits a continuous 
family of SAE. This is a limit case in the sense that beyond this range the operator 
becomes unbounded below. It has the peculiarity that only one SAE is locally scale 
invariant near the singularity at x — 0. This is reflected in even more pathological 
properties of its associated spectral functions. In Section H] we briefly review the 
singularities of the (- function, which also presents a branch cut, and the behavior 
of the heat-kernel trace, which in the general case admits a small-t asymptotic 
expansion in terms of negative integer powers of log t. 

In Section [5] we consider a locally homogeneous second order differential operator 
on the half- line R + [35]. In order to get discrete spectra for the SAE of this oper- 
ator we also include a quadratic term in the potential, V(x) = x 2 + [y 2 — |) jx 2 . 
We obtain similar results as in the compact case discussed in Section |3j Only two 
of the continuous family of SAE admitted by this Schrodinger operator for certain 
range of the coupling [5T] have domains which remain invariant under scaling trans- 
formations. The associated spectral functions show the same properties as for the 
regular potential case. On the contrary, the conditions satisfied near the singularity 
by the functions belonging to the domains of definition of the remaining SAE intro- 
duce a dimensionful parameter which explicitly break this local scale homogeneity, 
which implies the unusual behavior of the (^-function poles and powers of t in the 
asymptotic expansion of the heat-kernel trace already described. These results are 
obtained employing the von Neumann's theory of selfadjoint extensions of symmet- 
ric operators [H] to characterize the spectrum of the general SAE. They are also 
confirmed by an argument based on the asymptotic growth of the eigenvalues. 

Let us mention that this model corresponds to a classically integrable system 
pf5] whose quantum spectrum -subject to Dirichlet boundary conditions at the 



: Sclfadjoint extensions for more general singular potentials have been studied in [521 155] . 
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singularity- has been determined in [19]. This operator also appears as an effec- 
tive radial Hamiltonian for an isotropic harmonic oscillator in multi-dimensional 
Euclidean space with given angular momentum. For this model several results are 
known: the heat-kernel and the resolvent have been determined (also imposing 
Dirichlet boundary conditions) in |l6l U7J HB] ; the resolvent for a different boundary 
condition can be obtained from the Dirichlet-case by means of the so called Krein's 
formula [l9j [50] , which will be examined in detail in Section [6j 

In SectionOwe derive a generalization of Krein's formula [51] which can be applied 
to the kind of locally homogeneous operators with singular coefficients considered 
in this review. Krein's formula establishes a relation between the resolvents corre- 
sponding to different selfadjoint extensions. We generalize this relation to the case 
of a Schrodinger operator with a potential that contains terms with the same scaling 
dimension as the derivative term. Our generalized Krein's formula will make man- 
ifest the presence of non-standard powers of A in the large- |A| expansion of these 
resolvents. 

2. The first-order operator and its selfadjoint extensions 

Following [36], in this section we consider first-order symmetric differential oper- 
ator 

(2.D D =(°Av)' 
where 

(2.2) A = -d x + — , A = d x + - 

x x 

and \g\ < 1/2, defined on the domain T>(D) := C 2 x Qj°(0, 1), i.e. the set of smooth 
two component functions with compact support in the open segment (0, 1). For the 
given range of the coupling g, D is not essentially selfadjoint, admitting a continuous 
family of SAE. 

We will show that, for a general SAE of D, the large-|A| asymptotic expansion 
of the resolvent-trace Ty{(D — A) -1 } presents non-standard powers of A whose ex- 
ponents depend on the parameter g. Besides, the corresponding (- and ^-functions 
present simple poles lying at (/-dependent points in the complex plane, with residues 
which depend on the considered SAE. 

We proceed as follows: firstly, we construct the resolvents for two particular SAE 
for which the behavior of the functions in the definition domain of the SAE near 
the singular point x = is invariant under the scaling x — > ex. The asymptotic 
expansion of the resolvent for these special extensions displays the usual powers of A, 
then leading to the usual poles for the corresponding (- and 77-functions. Secondly, 
we show that the resolvents for the remaining SAE are convex linear combinations of 
these special extensions with A-dependent coefficients. This additional dependence 
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on A leads to non-standard powers in the asymptotic expansion of the resolvent, and 
hence to non-standard poles for the (- and ^-functions. 

These SAE are not locally invariant under the scaling x — > cx in the sense that 
the conditions defining the behavior of the functions belonging to its domain near 
x = are not scale- invariant. Rather, as c — > they tend to the conditions for the 
domain of one of the locally scale invariant SAE, while as c — » oo they tend to the 
other. As a consequence, the residues at the anomalous poles of the ^-function tend 
to zero as c — > and diverge as c — > oo. This behavior will be explained by means 
of a scaling argument. 

Integration by parts shows that the operator D is symmetric on T>(D). The 
adjoint operator D\ which is the maximal extension of D, is defined on the domain 

V(D^) of functions = ( J ^ ^ X ^(O) having a locally sumable first 

derivative and such that 

Since D is symmetric, T>(D) C T>(D^). The following lemma characterizes the 
behavior at the singular point x = of the two components of the functions in 
£>(£>t). 

Lemma 2.1. If &(x) E V(D^) and —~<g<\, then 
(2.4) | (t> x (x) - d[$] x 9 \ + | <p 2 {x) - C 2 [$] x~ 9 \ < K g \\D${x)\\ x 1 ' 2 , 
for some constants K g , Ci[$] and C 2 ^&\, where || ■ || is the L2 -norm. 
Proof: Eq. (12.31) implies 



Jo 



(2.5) 0x(x) 

(2.6) <t> 2 (x) =C 2 [<S>}x- a + x- 9 I y 9 fi(y)dy 
and taking into account that 



n 



(2.7) 



{21 



y 9 h{y)dy 



y 9 h{y)dy 



x <?+l/2 



<-=_||A| 

V 1 + 2 9 

r -S+l/2 

< ^7f=^= II/2 
V 1 - 2 9 



we immediately get Eq. (TjOJ with K g = (1 - 2a)~ 1/2 + (1 + 2a)~ 1/2 . 

The following lemma will be useful to describe the selfadjoint extensions of D. 



□ 
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Lemma 2.2. Let $(x) = ( ) ,*(x) = ( ^A X \ ] 



eV(D*). Then 




= {d[*]*C 2 [$] - C 2 [*]*d[$]} + {^ 2 (l)*0i(l) - Vi(l)*0 2 (l)} 




from which, taking into account the results in Lemma |2~TI Eq. (12.91) follows directly. 



Next, if \I/(x) in Eq. (12.91) belongs to the domain of the closure of D, D = (£)'), 
i. e. 



then the right hand side of Eq. (12. 9 1) must vanish for any <f>(x) € D(D^). Therefore 



On the other hand, if \l/(x) and <&(x) belong to the domain of a symmetric extension 
of D -which must be contained in T>(D*)- then the right hand side of Eq. (12. 9p must 
also vanish. Thus, each closed extension of D corresponds to a subspace of C 4 under 
the map $ — > (Ci[$], C 2 [$], 0i(l), 2 (1))- If we define the orthogonal complement in 
terms of the symplectic form of the right hand side of Eq. (I2.9p . then the selfadjoint 
extensions correspond to those subspaces ScC 4 such that S = S- 1 . 

Since we are interested in the consequences of the singularity at the origin, for 
simplicity in the following we will consider SAE satisfying the local boundary con- 
dition 



□ 



(2.12) 



*(ar) E V{D) C Dpt) , 



(2.13) 



C x \9\ = C 2 [tf] = tf(l) = 0. 



(2.14) 



0i(l) = O. 



Each extension is then determined by a condition of the form 



(2.15) 



aCi[$] +(3C 2 [$] = 



with a, (3 G K and a 2 + (3 2 = 1. We denote this selfadjoint extension by D^ a '^\ 
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2.1. The spectrum. In order to determine the spectrum of an arbitrary selfadjoint 
extension D^- a,l3 \ we need the formal solutions of 

( A(j> 2 {x) = \<f>i{x) , 
(2.16) (D - A)$(x) = =>- < 

[ A<pi(x) = X(f)2(x) , 

satisfying the boundary conditions given by eqs. ( 12. 14ft and ( 12.151) . 
The zero-mode, corresponding to A = 0, is given by 

dx 9 



(2.17) $(x) 



CoX-3 



The boundary condition at x = 1 chosen by Eq. (12. 14|) implies C\ = 0. The boundary 
condition at x = given by Eq. (12 . 1 5[) tells that C 2 = unless (3 = 0. Consequently, 
only the selfadjoint extension D^' ^ has a zero-mode. 

For A ^ we apply A to the second line in Eq. ( 12.161) so that, using the first line 
in Eq. (12.161) . one gets 

(2.18) {fl£_£^ + = 

The solutions of this equation take the form 

(2.19) 0i(x) = K 1 y/XJi_ g (X) + K 2 y/XJ g _i(X), 

where we have made the rescaling X := Xx and defined A = +vA^; K\,K 2 are 
complex constants. The lower component results 

(2.20) fa(x) = a {-^VX J_ g _i(X) + K 2 VX J g+ i(X)} , 
where a = A/A. 

From the behavior of <f>i(x) , 4> 2 (x) at the origin we can get Ci [$] , C* 2 [$] , respec- 
tively. The boundary condition now reads 

a K 2 X 9 3 Ki X~ 9 

(2.21) a C^] + 8 CM = , 2 ~ V g -iWi ^ = ° • 

2 9 2T (i + g) 2 9 2T (i -g) 

Let us first consider a particular selfadjoint extension, namely a = and 3 = 1. 
For a = 0, Eq. (12.21)) implies K\ = 0. Therefore, from the condition 0i(l) = we 
obtain the spectrum equation J g _i(X) = 0. Thus the spectrum of D^' 1 ^ is given by 

(2.22) A±, n = ±jW n = l,2,... 

y 2 

where is the n-th positive zero of the Bessel function J u (z). The spectrum is 
non-degenerate and symmetric with respect to the origin. 
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Figure 1. Plot for F(X) := X 29 



for g = 1/3 and p(a, (3) = 3. 



For a^O, Eq. ( I2.2ip can be written as 

4 9 r(| + s) 



(2.23) 



k; = cxX 



r (§-<?) 



/3 



In this case, the boundary condition at x = 1 determines the eigenvalues as the 
solutions of the transcendental equation 

X*-LJ± s L = ap (a,p), 



(2.24) 

where we have defined 

(2-25) p(a,/3):- ,. , , , 

For the positive eigenvalues A = A, o = 1 and Eq. (12.241) reduces to 



4 9 r(| + s) //3 



(2.26) 



F(X):=X 2 *-pj^ = p(a,(3). 

J Q-±{X) 



This expression has been plotted in Figured] for a particular value of p(a, 0) and g. 

On the other hand, for negative eigenvalues A = e ln X, a = e~ l7T and Eq. (I2.24p 
reads 



(2.27) 



F(A) = e- i >(ft,/3) = p(ft,-/3). 
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Therefore, the negative eigenvalues of D^ a '^ have the same absolute value as the 
positive eigenvalues of D^'^. 

Notice that for any selfadjoint extension the spectrum is non-degenerate and that 

there is a positive eigenvalue between each pair of consecutive zeroes of J„_i(A). 

y 2 

Moreover, the spectrum is symmetric with respect to the origin only for the selfad- 
joint extension corresponding to a = (from now on, the "D-extension" , see Eq. 
( I2.22p ) and for the extension corresponding to = (which we will call the "N- 
extension"). Indeed, from eqs. ( I2.24p and ( I2.25|) one can see that the eigenvalues of 
D^' ^ are given by 

(2.28) A = 0, \ ± , n = ±j ( rL, n = l,2,... 

2 y 

2.2. The resolvent. In this section we will construct the resolvent of D 

(2.29) G{\) = {D - X)- 1 , 

for its different selfadjoint extensions. 

We will first consider the two limiting cases in Eq. (I2.15p . namely the "D-extension" 
and the ' 'iV-extension" . The resolvent for a general selfadjoint extension will be later 
evaluated as a linear combination of those obtained for these two limiting cases. 

The kernel of the resolvent 

(2.30) G(x,y;X)=( ^f>^| ^f^^) 
K J V ; V G 2 i(x,y;X) G 22 (x,y;X) J 

satisfies 

(2.31) (D-X)G(x,y;X) = 6(x,y)l 2 , 
from which we straightforwardly get for the diagonal elements 

\dl - + A 2 1 G n (x, y- X) = -X 5(x, y) , 

(2.32) 

\dl _ ^tll + A 2| G22 ( X; y . X) = -x 6(x, y) , 

while for the non diagonal ones we have 

G 21 (x, y; A) = - {-d x + ^} G u (x, y; A) , 

(2.33) 

G 12 (x,y; A) = - [d x + ^| G 22 (x,y; A) , 

assuming A ^ 0. 



SPECTRAL FUNCTIONS OF NON ESSENTIALLY SELFADJOINT OPERATORS 



13 



Since the resolvent is analytic in A, it is sufficient to evaluate it on the open right 
half-plane. In order to do that, we use the upper and lower components of some 
particular solutions of the homogeneous equation (12.161) . 

Let us define 



(2.34) 



' L?(X) = 


Vxj g _ h (x), 


L°(X) = 


Vxj g+1 _(x), 


L» (X) = 


VxJ^ g (x), 


L»{X) = 


Vxj_ g _ h (x), 



R X {X- A) = VX J x (A) JiJX) - Jx (X) J _ 4 (X) 



R 2 (X; A) = y/X J g _i(A)J „ x(X) + Jx (A) J. +h (X) 



Notice that i?i(A;A) = AR 2 (\x;\) 



0. We will also need the Wronskians 



x=l 



(2.35) 



W [L^(X),R 1 (X;X)] = — cos(<?7r) J ff _x(A) 
W [L?(X), R 2 (X- A)] = - cos^ tt) J g _x (A) 

L 7T 



W [L?(X),i2i(X;A)] = -- cos(^Tr) Jx_ s (A) 

W [L%{X),R 2 (X;X)] = -~ cos(<?7r) Jx_ fl (A) 
which vanish only at the zeroes of J^(A), for v = ± (| — g). 



7d(A) 
-1 

1 

' 7iv(A) 
1 

' 7iv(A) 



2.3. The resolvent for the /J-extension. In this case, the function 
(2.36) = J*G D (x,y;\) ( f f f^ ) dy 

must satisfy the boundary conditions 0i(l) = and C 2 [$] = 0, for any functions 
A(x),/ 2 (x)GL 2 (0,l). 
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This requires that 

Lf(X)R 1 (Y;X), forx<y, 



(2.37) G°(x,y;X)= lD (X)x 
and 

(2.38) Gg(x } y;\) = - lD (\)x 



R 1 (X;X)L?(Y), ioix>y, 

L®{X) R 2 (Y; A), for x < y. 
R 2 {X;X)L°{Y), ioix>y. 



whereas the components G X2 (x, y\ A) and G 21 (x,y; A) are given by Eq. f |2.33|) . 
The fact that the boundary conditions are satisfied, as well as (D — A) = 
f x (x)\ 

f (x)J ' Can ^ e s t ra ishtforwardly verified from eqs. (I2.34p and ( I2.35p . Indeed, from 
eqs. qZM-EM one gets 

(2.39) faix) = C? [$]x 9 + 0(v^), <h(x) = O(V^) , 

with 

-ttA^ 1 



Cf[d>] 



2l+«cos(^7r)J g _i(A)r(| + ^) 
(2.40) 



x 







R l (Xy;X)f 1 (y)-R 2 (Xy;X)f 2 (y) 



dy, 



for A not a zero of J n _i(X). 

y 2 

Notice that Cf[$] ^ if the integral in the right hand side of Eq. (I2.40p is non 
vanishing. 

2.4. The resolvent for the iV-extension. In this case, the function 
(2.41) = jT G N (x, y; A) ( ) rfy 

satisfies the boundary conditions 0i(l) = and Ci[$] = 0, for any functions 
h(x),f 2 (x) GL 2 (0,1). 
This requires that 



(2.42) G? 1 (x,y;X)= lN (X)x 



Lf (X) i?i(F; A), for x < 
i?i(X; A)Lf (F), for x > y. 
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and 

L 2 y (X)R 2 (Y; A), for x < y, 



2.43) Gg(x,y;X)= lN (X)x 



R 2 (X;X)L%(Y), forx>y. 



whereas the components G^ 2 (x, y; A) and G^x, y\ A) are given by Eq. ( 12.330 . 

ffJx) 

These boundary conditions, as well as the fact that (D — A) = ( j )l 

can be straightforwardly verified from eqs. ( 12.341) and (12. 35 p . In this case, from eqs 
(I2T331 - f2735|) and fl2^D- l2^3j) one gets 

(2.44) <h(x) = 0{y/x) , 4> 2 {x) = C?l$]x- 9 + 0(V^), 

with 



(2.45) 



c?m = ^— ^ x 

2*-'c B(gir)Ji_ g {\)r(±-g) 



x 



Ri(\y,\)A(y) - R2(\y,X)f2(y) 



dy 



for A not a zero of JiJX). 

Notice that C 2 [®] ^ if the integral in the right hand side of Eq. ( 12.450 (the 
same integral as the one appearing in the D-extension, Eq. ( 12.400 ) is non vanishing. 

2.5. The resolvent for a general selfadjoint extension of D^ a ^ . For the gen- 
eral case, we can implement the boundary conditions 

(2.46) <M1) = 0, aCi[$] + /3C 2 [$] = 0, a, (3^0 



on 



h(y) 



(2.47) = J G(x,y;X)[ J f ^ j dy 

for any fi(x),f 2 (x) G L 2 (0,1) by taking a linear combination of the D- and N- 
resolvents 

(2.48) G(x, y; A) = [1 - r(A)] G D (x, y; A) + r(A) G N (x, y; A) . 

Since the boundary condition at x — 1 is already satisfied, we must now impose 

(2.49) a [1 - r(A)] C?[$] + r(A) C»[Q] = . 
Notice that, by virtue of eqs. ( gZQD , (12^451) and (l2~26l) . 

(2.50) aCf[$] -/3Cf[$] = 
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whenever A is an eigenvalue of D^ a,l3 \ Therefore, from Eq. (I2.49p we get the resolvent 
of D( a ^ by setting 

r(A)- ~ 1 - 



aCf[$]-/3Cf[$] p(a,/3) 

F(A) 

(2.51) 



1 



1 X 29 J ±- 9 ( x ) 



for A not a zero of J„_ i(A). 

y 2 

2.6. The trace of the resolvent. It follows from Eq. (I2.48[) that the resolvent of 
a general selfadjoint extension of D can be expressed in terms of the resolvents of 
the two limiting cases, Go (A) and Gn(X). Moreover, since the eigenvalues of any 
extension grow linearly with n (see section I27T1) . these resolvents are Hilbert-Schmidt 
operators and their A-derivatives are trace class. 
From the relation 

G(A) 2 = d x G(X)=d x G D (X)- 

(2.52) 

-r'(A) [G D (X) - G N (X)} - r(A) [d x G D (X) - d x G N (X)} , 

it follows that the difference Go (A) — Gn{X) is a strongly analytic function of A 
(except at the zeroes of t'(A)) taking values in the trace class operators ideal. 

From the explicit expressions of G D (X) and G N (X) (see Eqs. f[2T3?j) . fT2T38|) . ( j232D 
and (I2.43P ) we straightforwardly get 



Tr{<9 A G D (A)} = f tr{d x G D {x,x;X)}dx 
Jo 



2 



Similarly, 

(2.54, ^ GdW _ CwW} = _^_^ + _^. 



Moreover, since 



(2.55) «9 A Tr{G D (A) - G N (X)} = Tr{d x G D (X) - d x G N (X)} , 
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we get 



(2.56) 



A 



2 A JiJX) 



Tr{d x G D (\)-d x G N (\)} 
2A + J i(A) 



Finally, we can also write 



(2.57) 



Tr{G(A) 2 } = Tr{9 A C7 D (A)} - d A r(A) Tr{C7 D (A) - C7 A r(A)} 



2.7. Asymptotic expansion for the trace of the resolvent. Using the Hankel 
asymptotic expansion for the Bessel functions [55] (see Appendix |B]), we get for the 
first term in the right hand side of Eq. (I2.57P 



TrRG D (A)}~£^ 



(2.58) 



9_ 

A 2 



i a 



k=2 

g(g-l) Sg(g-l) 



+i a 



A 3 2 

(fl-3) (g-l)g(g + 2) 
2 A 5 



A 4 



+ol - 



where a = 1 for 3(A) > 0, and a = —1 for 9?(A) < 0. The coefficients in this 
series can be straightforwardly evaluated from eqs. (I2.53|) and (lB.14j) . Notice that 
A k (g, -1) = A k (g, 1)*, since A 2 k{g 1 1 ) is r eal an d A 2k +i(g, 1) is purely imaginary. 
Similarly, from Eqs. fl2~34l) . ( 1236]) and (lB~T5l) we get 

(2.59) Tr{G D (\)-G N (\)}~~ 9 



(2.60) 



Tr{d x G D (\) - d x G N {\)} ~ JL . 



A 

A 2 



On the other hand, taking into account Eq. flB.7j) . we have 

k 



(2.61) 



r(A) ~ < 



E 

k=l 



for 



1 

- < q < 
2 y 



p(a,/3) 

^(p(«,/3)e— (|-«)A- 2ff ) , for0<^<- 



fc=0 



where a = 1 (er = —1) corresponds to 3(A) > (3(A) < 0). Notice the appearance 
of non integer, ^-dependent powers of A in this asymptotic expansion. 
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Similarly 



(2.62) r'(A) ~ < 



2g , / e ain ^~ 9) X 2g 



fc=l 



p{a,0) 



for < o < 

2 y 



9 00 fc 1 

~T ^ * ( P( "' ^ ^"^^ ' , for < ^ < - , 

fc=i 



which are the term by term derivatives of the corresponding asymptotic series in 
Eq. (JZ6U). 

Collecting these results, we have 



Ox 



r(X)Tr{G D (X)-G N (X)} 



(2.63) 



SZ. /V^(§-s) \ I 

'2gk- l)X 29k - 2 , for - - < g < 



2gJ2 (p{a^)e-^^) (2 5 fc + 1) A" 2 '*" 2 , for < 5 < - . 



k=0 



Notice the ^-dependent powers of A appearing in these asymptotic expansions. 



2.8. The ( and 77 functions. The ^-function for a general selfadjoint extension of 
D can be defined, for 3?(s) > 1, as [9] 



(2.64) 



CM 



1 /" A 



2iti J c s 



-Tr {G(X) 2 } dX 



where the curve C encircles counterclockwise the spectrum of D. According to Eq. 
( 12.57!) . we have 



(2.65) C{s )=C D (s) + 



X 



l-s 



27T i Tr s — 1 



r(X)Tr{G D (X)-G N (X)} 



dX 



where ( D {s) is the (^-function for the D-extension. 

Since the negative eigenvalues of the selfadjoint extension D^*'^ are minus the 
positive eigenvalues corresponding to extension _d£ q ' - ^ (as discussed in Section [27Tj) , 
we define a partial (^-function by means of a path of integration that encircles the 



SPECTRAL FUNCTIONS OF NON ESSENTIALLY SELFADJOINT OPERATORS 19 

positive eigenvalues only, 

i pioo+O \1— s 

(2.66) 

rioo+O r- 

<7A 



= Cf 00 - / ^— r 5 A r(A) Tr{G D (A) - G^A)} 
where C+ ( s ) is the partial (^-function for the D-extension. Notice that we can write 

(2.67) 

where h(s) is an entire function. In order to determine the poles of £+ (s), we add 
and subtract a partial sum of the asymptotic expansion of Tr |G(A) 2 } (obtained in 
Section [2~7|) to the integrands of the right hand side of Eq. (12.671) . 
In this way, we get for the D-extension and for real s > 1 



A' 



K ' k=0 

where h^(s) is an analytic function in the open half plane 9ft(s) > 1 — N. In 
consequence, the meromorphic extension of C+ ( s ) presents a simple pole at s = 1 
(see Eq. (12.681) ) with a residue given by (see Eq. (12.531) ) 

1 r°° +0 f ^„+i(a)] 1 

(, 69) ResCfWU ._/_ +/ o aA (^Uj, A _. 

It also presents simple poles at s — —k, for k — 0,1,2, ... , with residues given by 



(2.70) ResC+(s 



%{A k+2 (g,l)} 



I s=— k 



7T 



(* + l) 



with the coefficients ^(p, 1) given by Eq. (12.581) . In particular, notice that these 
residues vanish for even k. 
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For a general selfadjoint extension Dx^\ we must also consider the singularities 
coming from the asymptotic expansion of d\[r(\) Tv{Gd(^) — Gn(\)}] (see Eq. 
(I2.63P ). For definiteness, let us consider in the following the case — ~ < g < (the 
case < g < | leads to similar results). 

From Eq. (12.661) . and taking into account Eq. (I2.67p . for real s > 1 we write 



N 

E 



9 



7T (s 



Hm(s) 



(2.71) 



tt(s-I) A 



e _ if( _ s _ 1)/xl _ s 



AT 

E 



e 2 



(2gk-l)fx 29k - 2 \ dfx 



29 



A' 



7T S 



i) E 

y fc=i 



2gk-l 
s — 2 g k 



5? 



if (fc-s-l) 



e 2 



p(«,/3) fe 



^TV(S) 



1 ' 



(«,/3), 



S - 



where Hn(s) is holomorphic for 3fJ(s) > 2g(N + 1). We conclude that £. 
C+ (s) has a meromorphic extension which presents a simple pole at s = 1, with a 
vanishing residue, 



8=1 



(2.72) 



2vri 



i oo+O 



Res (C^\s) - C?00 

V(A)Tr{G D (A) -Gat(A)} dA = , 



A°<9> 



i oo+O 



as follows from Eqs. (1239|) and (l2lJTj) . 

Notice also the presence of simple poles located at negative non integer ^-dependent 
positions s = 2gk for k = 1, 2, . . . , with residues which also depend on the selfadjoint 
extension, given by 



(2.73) Res{#^(*)-C?(a)} 



-2s 



sin 



- - ) A- 7T 



Following the discussion after Eq. f)2.27p . we get for the complete ^-function 



(2.74) 



C 



(«,/3), 



c 



c 



In particular, for the D-extension -whose spectrum is symmetric with respect to the 
origin (see Eq. f 1 2 . 2 2 [) ) - we get 

(2.75) C D (s) = (l + e- i7rs )C?(s). 
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Notice that ( D (s) is an entire function. Indeed, from Eq. ( 12.70ft . the residue at 
s = —k vanishes for k even, whereas for k — 2 1 + 1, with Z = 0, 1, 2, . . . , we get 



(2.76) 



Res {C D 0)}| 



-21-1 



(l + e"( 2!+1 ») Res {(°(s)} 



s=l-2l 



0. 



On the other hand, for a general selfadjoint extension, the sing ularities of ( {a ^{s) 
are simple poles located at s — 2 g k < 0, for k — 1, 2, . . . , with residues 



Res {C 



C D (s)} 



s=2 g k 



(2.77) 



Res 



+ e~ 



C?00 



_^kll sin ( 2 9 h*) ci7r (i- g ) k 
' 7T p(a,/3) k 



where we have used p(a, —0) = —p(a,(3) (see Eq. (12.251) ). 



s=2 g k 



Similarly, for the spectral asymmetry [56J we have 



(2.78) 



(«-/9), 



In particular, ^'^(s) = r]^ 1 '°- ) (s) = 0, since the corresponding spectra are symmetric 
(see eqs. (I2.22f) and (12.281) ). For a general selfadjoint extension and — | < g < the 
function r)( a '^(s) presents simple poles at s = 2gk, for odd k = 1,3,5,..., with 
residues given by 

Ag_ sin [(\-g) kn] 

71 



(2.79) 



Res {v (a 'P\s)} 



s=2 g k 



P (a,py 



For the case < g < |, an entirely similar calculation shows that C+ (s) — C+ ( s ) 
has a meromorphic extension which presents simple poles at negative non integer 



^-dependent positions, s = — 2 g k, for k = 1,2, 
selfadjoint extension, given by 



Res {C 



(«,/3), 



c 



Di 



with residues depending on the 



-2gk 



(2.80) 



- — p(a,(3) k sin 

71 



9 kn 



From this result, it is immediate to get the residues for the (- and 77-functions. One 
gets the same expressions as in the right hand sides of eqs. (12. 77ft and (12.791) . with 

p(a,/3) and e ln ^~ 9 ^ k replaced by their inverses. 
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2.9. Scale Invariance. Let us remark that when neither a nor f3 is 0, the residue 
of C+ at s = —2 \g\ A; is a constant times (/3/a) fcsign< - 9 - ) . This is consistent with the 
behavior of D under the scaling isometry Tl u(x) := L^ 1 / 2 u(x/L) mapping L2(0, 1) 
onto L 2 (0,L). The extension D^ a '^ is unitarily equivalent to the operator L D^ a '^ 
defined on L 2 (0, L), with a := L g a and $ := L~ g f3: 

(2.81) T L D {a ' p) := L D^' $) T L . 

Notice that only for the extensions with a = or (3 = the boundary condition at 
the singular point x = -given by Eq. f !2. 1 5[) - is scale invariant. 
The partial (^-function of the scaled operator can be written as 

(2.82) (^\s) = L s (^\s), 
whereas its residues satisfy 

(2.83) Res (cf 4) (s)) = L~ 2|9|fc Res 

P + ) s=-2\g\k I J s=-2| S |fe 

The factor L _2 ' 9 ' fc exactly cancels the effect the change in the boundary condition 
at the singularity has on p(a, (3), 

(2.84) p(a, p) ksi ^) = L 2 lsl k p{a, p) ksi ^) _ 

Thus the length of the intervals (0, 1) and (0, L) has no effect on the structure of 
these residues, which presumably are determined locally in a neighborhood of x = 0. 

Finally, let us point out that these anomalous poles are not present in the g = 
case. Indeed, in this case r(A) in Eq. (I2.5ip has a constant asymptotic expansion, 
whereas Tr{Go(A) — Gat(A)} ~ (see Eq. (12.591) ) . Moreover, the residues of the 
poles coming from C+ ( s ) vanish (see Eqs. ( 12.70P and (12. 58ft ). except for the one at 
s — 1, whose residue is 1/ir (see Eq. (I2.69P ). 

Consequently, the presence of poles in the spectral functions located at non integer 
positions is a consequence of the singular behavior of the 0-th order term in D near 
the origin, together with a boundary condition which breaks the scale invariance. 

3. The second order case 

In this section we will show that one obtains similar results for the selfadjoint 
extensions of the second order differential operator 

(3.1) H = -g i+ ^, 

x z 

defined on a set of functions (f>(x) G C^°(0, 1). Since we are interested in the effects 
of the singularity at x — 0, for definiteness we impose again 0(1) = 0. Proceeding 
as in the previous section [37], for \g\ < | (to be congruent with the first-orden 
case analyzed in Section [2]), one can show that H admits a continuous family of 
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selfadjoint extensions in L 2 (0, 1), H^ a, ^\ characterized by two real parameters a, (3 
satisfying a 2 + (3 2 = 1 and defined on functions which behave near the origin as 

(3.2) <p(x) = C 1 x 9 + C 2 x 1 - 9 + 0(x i/2 ), 
where the coefficients Ci,C 2 satisfy 

(3.3) ad + f3C 2 = 0. 

The spectrum of H^ a '^ is determined by the relation (analogous to Eq. ( 12. 26ft ) 

(3.4) ^(fi):=-F( f j l ) = Q(a,P), 
where now we define 

(3,) e(tt ,^g) 2 ^E|±|. 

Also in this case, a = and = correspond to two scale invariant boundary 
conditions at the singularity. For these two limiting extensions, it is easily seen from 
Eqs. ( ESI) , A2~39l) and that the resolvent of satisfies 

(3.6) g D ' N (x,y- f i 2 ) = -G^ N (x,y ]f i). 

The resolvent for a general selfadjoint extension H^ a '^ is constructed as a convex 
linear combination of Q D (fi 2 ) and G N (fi 2 ) as in ( 12.481) . with a coefficient 

(3-7) r(/i) = ! R , ■ 

Following the methods employed for the first order case [37] , one can show that the 
^-function associated to H^ a '^ has a meromorphic extension which presents simple 
poles located at negative (^-dependent positions, 

(3.8) s k = -(^-g)k, for fc = l,2,... , 

with residues which depend on the SAE given by 

Res {((^)( s )-C D (s)}\ s=Sk = 

(3.9) 

' i)0( a ,/3)*sin '-{2g-l)k 

Notice that these poles are irrational for irrational values of g. Moreover, these 
residues vanish for the "N-extension" (g(a,0) = 0), and have a singular limit for 
a — » 0. As in the first order case, this behavior can also be explained through scaling 
arguments 
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Finally, let us remark that the relation between the (^-function and the trace of 
the heat-kernel of D^*'^ (See Appendix [A} straightforwardly lead to the following 
small-t asymptotic expansion, 



Tr 



-tD 



(«,/3) 



-tDl 



1 

9 ~2 



(3.10) 



k=l 



k 



\^ Q {a^) k sin \l (2g - 1) k 

A 71 12 



-g)k 



The first term in the right hand side coincides with the result reported in [34J . Notice 
also the (^-dependent powers of t appearing in the asymptotic series in the right hand 
side of Eq. (I3.10P for the general SAE and the dependence of the SDW coefficients 
on the (a, /3)-parameters. In particular, the first term in this series reduces to 

(3.1!) 



a Y{\ 



9) 



This power of t also coincides with the result quoted in 
coefficient. 



but we find a different 



4. The g = ± case 
The second order differential operator given by Eq. (13. ip with g 

d 2 1 



i 

2' 



(4 ' 1} Ho = 'd^ 4x2' 

defined on the domain <f>(x) G C^°(0, 1), also admits a continuous family of nontriv- 
ial SAE in L2(0, 1), which reflects in even more unusual properties of its spectral 
functions. This problem has been considered by K. Kirsten et al. in |40 | |4T | I42 | |4"3], 
articles where an error in Appendix A of [2Z] has been corrected. 

Also here, for definiteness, we impose (f)(1) = 0. As before, the SAE are 
characterized by two real parameters satisfying a 2 + /3 2 = 1 and defined on a domain 
of functions which behave near the origin as 

(4.2) <f>{x) = C iy ^ + C 2 ^\ogx + 0(x 3/2 ), 

where the coefficients Ci,C 2 satisfy Eq. (13.31) . 

The eigenfunction of := H^ ' 1 ^ corresponding to an eigenvalue A = fi 2 is given 
by, 

(4.3) (j)(x) = a/x Jq(llx) . 

The condition 0(1) = tells that jj, is a (positive) zero of the Bessel function Jq(z). 
On the other hand, for an arbitrary SAE H^*'^ with a 0, the eigenfunction 
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corresponding to an eigenvalue A = /i 2 is given by 

(4.4) <f>{x) = {d - C 2 (log/x/2 + 7)} J (»x) + I C 2 ^N (iix) , 

where Ci,C 2 are constrained by Eq. (13.31) . The condition (f)(1) = leads to the 
equation 

(4.5) 2(9 - log/i) JoGu) + 7riV (/x) = , 

where 6 1 := —/3/a + log 2 — 7, which determines the spectrum of Hq . Notice that 
there are no negative eigenvalues. 

The trace of the resolvent Gq([i 2 ) := (Hq — fj 2 )" 1 can be readily computed to get 

This trace admits the following asymptotic expansion in integer powers of /i 

zcr 1 icr 1 _ 5 

where a = +1 (—1) for > (S(/i) < 0). From this asymptotic expansion one 
concludes that the poles of the corresponding ^-function Tr {Hq} s are located at 
s = 1/2 -k, for k = 0,1,2,... 

On the other hand, the trace of the resolvent Gq* {l*?) '■= (-Hq — A 42 ) -1 ) corre- 
sponding to a general selfadjoint extension, gives 

(48) Tr i^° (/i) )-2^2(^-log At )Jo( Ai ) + vriV ( / ,) + 



^{kN^) + 2(9- log/i)] J (/i)' 



The asymptotic expansion of the first term in the R.H.S. of Eq. ( 14. 8 p is also given 
by expression (14.71) . Therefore, it leads to the standard poles of the corresponding 

f-function, Tr {Hry, located at .= !/*- ,, for * = 0, 1,2, . . . However, the 
asymptotic expansion of the second term in the R.H.S. of Eq. (14 ,8p is given by 

(4.9) JM - 1 



(j?[irN {n) + 2(9 - log/i) J Q {n)] y?\ma + 2(9 - log/i)] 
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with a = 1 (er = — 1) for /x in the upper (lower) half-plane. This expression gives an 
additional contribution to the ^-function given by 

Tr [h^Y* -Tr {H?}" = 
(4.10) =- — - (■'""" I ( 0. [ — 20 I s I — e _i7rs r (o, ( ZlJL — 26 



e- 2ds 

~ 2ni [' V" 7 V 2 

= _I e 2s (-- log2+7 ) sin(vrs) \ogs + H(s), 

7T 

where H(s) is an entire function of s. Therefore, the ^-function of the general SAE 
of H develops a cut on the negative real axis, which is present even for /3 = and 
has a singular behavior for a — > 0. 

As a consequence, the derivative of the (^-function in a neighborhood of the origin 
behaves as 



(4.11) A Tr {tf^} S = -l ogs + 0(l). 

This logarithmic branch cut from the origin for the general SAE makes the functional 
determinant in Eq. (\1.6\i ill-defined (See the discussions in jlQl Ej). Notice that 
this logarithmic cut is absent only for the unique SAE of Hq which is locally scale 
invariant, namely = H^ 0,1 \ 

On the other hand, the term (14.91) gives the following contribution to the heat- 
kernel trace of a general SAE of H 

Tr <^ e~ tH ° [ - Tr \ e~ tH ° 



(4.12) 



e 



iri/j, TcN (jj,) + 2(6 - log/i) J (aO 



1 r°° e~ tx 

- $s / — -dx + R(t) 

7T J l X (log X - 26 + Z7r) 



where i?(t) is a smooth function at t — 0. It was shown in [ID] that the integral in 
the right hand side of Eq. (I4.12p has an asymptotic expansion for small t in terms 
of negative integer powers of logt. 



5. Non-compact case 

In this Section we will consider a locally homogeneous (near a singularity) symmet- 
ric second order differential operator on the noncompact one-dimensional manifold 
M = M + . We will see that also in this case the associated ^-function presents a 
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non-standard singularity structure, related with the breaking of this scale homo- 
geneity by the definition domains of the SAE of this operator. We will employ 
a different approach to this problem, based on the von Neumann's theory of self 
adjoint extensions of symmetric operators. This allows to express the spectrum of 
the SAE in terms of a trascendental equation from which we are able to derive an 
asymptotic expansion of the eigenvalues. This eventually leads to the pole structure 
of the ^-function we are interested in. 

Then, following [22], let us consider the operator 



(5.1) 



with 



(5.2) 



H 



V{x) 



d 2 



dx 2 



V(x) 



+ x 2 , 



densely defined on the domain T>(H) = C^°(IR + ), the linear space of smooth functions 
4>(x) with x G R + and compact support out of the origin. We have added to the 
singular at the origin term of the potential an x 2 term in order to get discrete spectra. 

According to von Neumann's theory of deficiency indices [H], in order to get 
the SAE of H in L2 (K + ) we need to compute the adjoint and determine the 
deficiency subspaces K± := Kernel [W =F i). The domain of W is the subspace of 
square integrable functions having an absolutely continuous first derivative and such 
that 



(5.3) 



Hty(x) = -ip"(x) + V{x)ip(x) G L 2 (M + ) . 



Notice that no boundary condition is imposed at x — 0. To compute the deficiency 
indices n± : = dim )C± of H we must solve the eigenvalue problem 

(5.4) #V A = + V(x)Mx) = A^a, 

for <p x e V(W) and A G C with imaginary part 3(A) 7^ 0. Let us define the 
parameter 

(5.5) a:=l + Z y. 



For the case < v < 1, i.e. 1/2 < a < 3/2, for any A G C Eq. (15. 4p has a unique 
nontrivial square-integrable solutions given by 



(5.6) 



, ,'2a+ 1 — A 1 2 
>a( > ) = x e 2 U ( ; a + -; x 2 



where U is the confluent hypergeometric function as defined in [55J. Thus, the de- 
ficiency subspaces K,± are one-dimensional, n± = 1, and H admits a one-parameter 
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family of which are in a one-to-one correspondence with the isometries from 

JC + onto /C_. The deficiency subspaces JC + and /C_ are spanned by (f> + := (j)\ =i and 
0_ := 4>\=-i = 4>\, respectively Each isometry W 7 : fC + — > /C_ can be characterized 
by a parameter 7 £ [0, it) defined by 

(5.7) U 



-2*7, 



7S 



Each isometry is identified with a selfadjoint extension if 7 , a closed restriction of 
W to a linear subspace 

(5.8) V{H y ) C V(H^) = V(H) © /C+ © £_ , 

where if is the closure of if. The domain D(if 7 ) is defined as the set of functions 
that can be written as 

(5.9) = o + ^(0+ + e- 2i7 0-), 

with 0o ^ T^{H) and A a complex constant, and the action of the operator if 7 by 

(5.10) H, /( j) = H ] (f) = H ] (f) Q + 1A (<f) + - e~ 2 ^(f)_) . 

Let us consider, for simplicity, the repulsive case {y > |), that is 1 < a < 3/2. In 
Appendix [Cl we show that ^o^) = o(x a ) and O ( X ) = o{x a ~ 1 ). Therefore, from Eq. 
(15.91) and Eq. ( 15 . 61) with A = ±2 we get the following condition for <p{x) £ D(if 7 ) 
near the origiru, 

(5.12) 9, log^(x) = 1Z2 _ 2 E||Z^ « . + , 

x T(a - |) cos (7 - 72) 

where we have defined 71 = arg {T[(— 2a + 3 — «)/4]} and 72 = arg {r[(2a + 1 — «)/4]}. 

The boundary condition specified in Eq. ( 15.121) characterizes the domain of the 
selfadjoint extension if 7 . In order to determine its spectrum, we select from the set 



3 This is in accordance to Weyl's criterion [53] according to which, for continuous V(x), H is 
essentially selfadjoint if and only if it is in the limit point case, both at infinity and at the origin. 

In addition, if V(x) > M > 0, for x large enough, then H is in the limit point case at infinity. 
In consequence, in the present case H is essentially selfadjoint if and only if it is in the limit point 
case at zero. 

In particular, for positive V{x) {i.e. v 2 > j), if V(x) > jx~ 2 for x sufficiently close to zero 
then H is in the limit point case at the origin. On the contrary, if V(x) < (| — e)x~ 2 , for some 
e > 0, then H is in the limit circle case at zero. 

4 Notice that 

(5.11) 4>+{x) = x 1 -* ( , 1 , r ("~^ rr, + O (x 2 )) +x a ( „ r ,^~ a ) . + O (x 2 )) . 

Then, it is easy to see that there are two locally scale invariant SAE for which the functions in their 
domains behave near the origin as x a and x 1_Q respectively. For any other SAE, the simultaneous 
presence of both powers of x in the boundary condition near the singularity breaks this local scale 
homogeneity 
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of eigenfunctions of given in Eq. (15. 6 j) those which satisfy Eq. ( 15. 12ft . This leads 
to the following transcendental equation for the eigenvalues A [35] , 

(5-13) *^ Y = / 3( 7> «) ) 

r(i-« - %) 

where we have defined the constants 

2a + 1 
k := — — G [3/4, 1) , 

(5.14) 4 

/%, «) := cos (7 - 71)/ cos (7 - 72) . 

Eq. (15.131) determines the discrete spectrum of the self adjoint extension characterized 
by the parameter 7. From now on we will refer to the SAE as identifying it 

by the value of/3GRU{— 00} defined above. 

As expected, the spectrum of is bounded from below; however it presents a 
negative eigenvalue for those selfadjoint extensions characterized by (3 > T(k)/T(1 — 
k) (even though the potential V(x) > 0). Moreover, there is no common lower 
bound; instead, any negative real number is in the spectrum of some selfadjoint 
extension. 

For any value of v G [|,1), there are two particular selfadjoint extensions for 
which the spectrum can be easily worked out (see Eq. (15.131) ): 

• For (3 = the spectrum is given by 

(5.15) A n = 4(n + l-«), n = 0,1,2,... 

• For (3 = — 00 the spectrum is given by 

(5.16) A n = 4(n + K), n = 0,1,2,... 

For any other value of /3, the eigenvalues also grow linearly with n, 

(5.17) 4(n - 1 + k) < A n < 4(n + k) . 

5.1. Pole structure of the (-function. We will now study the pole structure of 
C^(s) corresponding to an arbitrary SAE H(py Notice that, since the eigenvalues 
grow linearly with n (see Eq. (15.171) ). (p^s) is analytic in the open half-plane 3?(s) > 
1. 

Let us begin by considering the (-functions for the SAE characterized by = 
and (3 = —00, which can be explicitly evaluated from the expression of its spectra 
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in Eqs. ( 15. 15ft and ( 15. 16ft . We obtain 

oo 

(5.18) C P=0 (s) = 4~ s £(n + 1 - = 4 "^(s> 1 " «) , 

?1=0 

OO 

(5.19) C^=-oo(s) = 4-* J> + K )~ s = 4~X(s, «) . 

n=0 

where C(s, (?) is the Hurwitz ^-function, whose analytic extension presents a unique 
simple pole at s — 1 with residue Res£(,s, q)\ s= i = 1. Therefore, for (3 = and 
/3 = — oo the (^-function presents a unique simple pole at s — 1, with residue 1/4. 
For finite /3, let us define the holomorphic function, 



recall that we are considering the repulsive case f < k < 1. The eigenvalues of the 



i 



selfadjoint operator correspond to the zeroes of /(A) which, consequently, are 
all real. They are also positive, with the only possible exception of the lowest one, 
according to the discussion in the previous section. Moreover, the zeroes of /(A) are 
simple and isolated; thus, the ^-function can be represented as the integral on the 
complex plane 



(5.2i) Ce ( s) = _L^ A -.ZW + e(-v„)A-j, 

where C is a curve which encloses counterclockwise the positive zeroes of /(A), A ,/3 
is the lowest eigenvalue and O(-) is the Heaviside function. 

For 5J(s) > 1 the path of integration in (" 15 .211) can be deformed to a vertical line 
to get 



<5 - 22) « S) = W^ + A ~W A + MSH 



where hi(s),h 2 (s) are holomorphic functions. Eq. (15 .22j) gives an integral represen- 
tation of C/?( s ); analytic in the half-plane 3?(s) > 1. To compute its meromorphic 
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extension to the whole complex s-plane and its pole structure we need the asymp- 
totic expansion of f'(X)/f(X), which is given by 

/'(A) 



1 l 00 

-log(-A) + Tc,( K )(-A)- fe + 



(5.23) 

oo oo 
JV=1 n=0 

The coefficients Cfc(ft) are polynomials in k. As we will see, these terms do not con- 
tribute to the pole structure of C,s( s )- On the other hand, the coefficients C^^k, (3) 
are defined through the following relations: 

(5.24) C N>n (K, (3) := - (4 2 - 1 /?) N (2k-1 + ^) b n {K, N) , 



N f 

00 ^00 

(5.25) N ) z ~ 2n '■= ex P { N E am ^ z ~ 2m 

n=0 ^ m=l 

(5.26) ^^^{[d-^-^^^iZ?) 



X 



[(1 - /t) 2m + /t 2m ] + (2m + 



X 



X ( 2p - 2 1 ) W [m " P)+1 ~ (1 - «:) 2(m - p)+1 ] } • 

Notice CW in (K, /?) = for /3 = 0. 

Replacing into Eq. (15.22)) the dominant logarithmic term in Eq. (15.23)) we get 
a simple pole at s = 1 with residue 1/4. The remaining terms in the asymptotic 
expansion of f'(X)/f(\) are of the form Aj(— A) - - 7 , for some j > (see Eq. (15.231) ). 
Replacing these terms into Eq. (15.221) gives simple poles at s — 1 — j, with residues 
given by — (Aj/tt) sin(7rj). Notice that these residues vanish for integer values of j. 

In conclusion, there is a simple pole at s = 1 with residue 1/4 as for the SAE with 
,5 = 0, —00. But for a general SAE with (3 7^ 0, —00 and for | < k < 1 there 
are also simple poles at non integer values of s. Indeed, for each pair of integers 
(N, n) with N = 1, 2, 3, . . . and n — 0, 1, 2, . . . the function C/?( s ) nas a simple pole 
at the negative value 

(5.27) sjv.n = -N(2k - 1) - 2n , 
with a /3-dependent residue given by 

(5.28) Res [&(*)] | s=Sjv _ n = C N , n { K , (3) sin(2vriV«:) . 
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Let us remark that when k is a rational number, there can be several (but a finite 
number of) pairs (N, n) contributing to the same pole. On the contrary, when k is 
irrational the poles coming from different pairs (N, n) -which are also irrational- do 
not coincide. 

Finally, notice that a pole of C/?( s ) & t a non integer sjv,n — —N(2n — 1) — 2n 
implies the presence of a term proportional to t N ( 2K - 1 )+ 2n j n the small-t asymptotic 
expansion of Ti{e- tH ^} (See Appendix |Aj) . 

Therefore, for this second order differential operator on the half-line we obtain 
similar results as in Section [3] for a compact segment: There are two SAE whose 
definition domains are locally scale invariant near the singularity (See footnote S]) 
and show the usual properties in their spectral functions. Moreover, there exists a 
continuous family of other SAE which presents anomalous poles in their (^-functions, 
dependent on an external parameter (the coupling u), with residues dependent on 
the SAE. 

We finish this section by presenting an alternative method based on the relation 
between the pole structure of the ^-function and the asymptotic growth of the 
eigenvalues, which confirms our results. Indeed, the pole structure of C^( s ) can a l so 
be obtained from the asymptotic expansion of the eigenvalues A n for n ^> 1. 

By solving Eq. (15.131) order by order in n we get [35J 

Ad 

(5.29) \ n = An + 4(1 - «) + — sin(2 tt k) n 1 ' 2 * + 

TT 

AB 2B 2 
+ — (1-3k + 2k 2 ) sin(2 7rfi;)^ 2K - — sin(4 tt k) n 2 ^ K + . . . , 

7T 7T 

where we have retained only powers of n greater than —2. This leads to 

(5.30) Cp(s) = 4" s C(s) + s A- s (k - 1) C(s + 1) + 

+s (s + l) 4" s ^-((s + 2)-s4- s - sin(2/t7r) ((s + 2k) - 

2 7T 

- s (s + 2 k) A~ s - (k - 1) sin(2 tt k) C(l + s + 2 k) + 
+s A~ s — sin(4 tt k) ((s - 1 + 4 k) + ... , 

2 71 

where ((z) is the Riemann (^-function. Taking into account that ((z) presents a 
unique simple pole at z = 1 with residue 1, Eq. ( 15.301) confirms, order by order in 
this development, the pole structure we have already found. 

6. Krein Formula 



In this Section we study the behavior of the resolvent of a locally homogeneous 
second-order differential operator in relation with its selfadjoint extensions, and the 
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consequences this has on the properties of other spectral functions. This will be 
done in the framework of the Krein's formula [54], which relates the resolvent of two 
SAE of the given operator. 

We consider the differential operator 

(6.1) A = -dl + ^^ + V(x), 

where v G (0, 1) C M and V(x) is an analytic function of x G K + . 

The operator (16. ip defined on T>(A) := C£°(M + ) admits a continuous family of 
selfadjoint extensions A 9 characterized by a real parameter which we call 9. Since 
the operator e~ tA corresponding to a general SAE A 6 is not trace class (notice that 
the base manifold M + is non-compact), we will consider the trace of the difference 
e~ tA<> — e~ tA °° , where A°° corresponds to the Friedrichs extension [33] . We will show 
in Theorem 16.151 that this trace admits an asymptotic expansion given by 

oo oo 

(6.2) Tr{e-" 9 -e-"~}~'Ea n { V ,V)ti + £ b N>n {v, V) 6 N ~\ . 

n=0 N,n=l 

As we will see, the SDW coefficients a n (u, V), 6jv,n(^, V) can be recursively computed 
for each given potential V(x). Let us remark that the singular term in (16.11) not only 
contributes to the coefficients a n (u,V) of the standard powers of t but also leads 
to the presence of non-standard powers of t whose exponents are not half-integers 
but depend on the "external" parameter v. We will also show that these terms are 
absent only for the SAE with 8 = and 9 = oo, which correspond to selfadjoint 
extensions characterized by scale invariant domains. In Section l6~4l we will consider 
the case V(x) = x 2 to compare with our results in Section [51 

Actually, the main content of this section is the derivation of a generalization of 
the Krein's formula to the case of these kind of operators with singular coefficients, 
from which the expansion (16.21) is obtained as a byproduct. We will make use of 
the two particular SAE of the operator A, namely A and A°°, for which the v 
dependent powers of t in (16.21) are absent. The expansion in Eq. (16.21) will come 
out as a consequence of the relation between the resolvents corresponding to an 
arbitrary SAE and to A°°. This relation is called Krein's formula [54] and has 
already been established for the case of operators with regular coefficients. We will 
therefore extend this result to A in (16. ip and then use this generalization to prove 
the asymptotic behavior (16.21) . 

6.1. The regular coefficients case. In this section we state without proof two 
theorems valid for the case of differential operators with regular coefficients. Theo- 
rem [6J] describes the selfadjoint extensions of a symmetric operator in terms of the 
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(regular) boundary values of the functions belonging to its domains [57J. The state- 
ment of theorem 16.31 is the Krein's formula (see [49]), which relates the resolvents 
corresponding to different selfadjoint extensions. 

Theorem 6.1. Let A be a symmetric operator densely defined on a subspace T>(A) 
of a Hilbert space H, for which the deficiency indices are equal n + = n_ =: n < oo. 
Then: 

• There exist two surjective maps T\, T 2 : V(A^) C n such that V <p, e 
P(At) 

(6.3) (0, Aty) H - (A*</>, ij) H = (r x 0, r 2 ^)c- - (r 2 0, r^k- , 

where (-, •)% is the inner product in 7-L and (-, -)c« is the usual inner product 
in C n . 

• The selfadjoint extensions A^ M > N ^ of A are characterized by two matrices 
M,N e C nxn , such that M ■ is hermitian and (M\N) e C nx2n has rank 
n. The domain of definition of A^ M,N > is defined as 

(6.4) V (A iM ' N) ) := {<j> e V(A ] ) : MT l( f) = NT 2 <j)} . 

Since the restrictions of I\, r 2 to Ker(A^ — A) are invertible, we can establish the 
following definitions: 

Definition 6.2. 

(6.5) r~ x (A) := (ri| Kc r(At-A)) 1 , 

(6.6) K(X) ^-^-^(A). 

Now we can write down the Krein's formula, which expresses the resolvent 
fj^(M,N) _ )\ corresponding to an arbitrary selfadjoint extension in terms of the 

resolvent (A^ 1 ' ^ — A) corresponding to the selfadjoint extension characterized by 
the matrices M — 1 and N = [58] . 

Theorem 6.3 (Krein's formula). 

(6.7) {A™ - A) -1 = (^.°> - A)" 1 + r r '(A) ■ (M + % m) ■ (^'(V)) 1 . 

Example: let us write down the Krein's formula for the case of the one-dimensional 
second order differential operator 

(6.8) A = -dl + U(x) , 

defined on Cg°(IR + ) C L 2 (R + ). We assume that the potential U(x) £ CQR+) is in the limit 
point case [33] at infinity and in the limit circle case [33] at x = 0. This is the case if there 
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exist 5, e > and a positive differentiable function f(x) > —U{x) such that 

(6.9) < U(x) < 3 ^ 4 ~ £ Vx G (0, S) , 

x z 

f°° dx fix) 

(6.10) / — diverges , and , is bounded near oo . 

Then, the deficiency indices of A are n± = 1 [44]. In relation to Theorem 16. 1} we define 
the boundary operators Li,r2 as 

(6.11) T x <t>{x) := 0(0) , 

(6.12) r 2 0(x) := 0'(O) . 

According to the second statement in Theorem 16. H the selfadjoint extensions A e of A are 
characterized by a real parameter 9 corresponding to M _1 iV £ 1 and their domains of 
definition V(A e ) are given by (see Eq. ([IT 



3.13) V (V) = {4> E V(A^) : 0'(O) - 0(0) = 0} . 



As expected, we obtain the classical boundary conditions of Robin type. The extension 
characterized by N = (Dirichlet boundary condition) corresponds to 9 = oo while M = 
(Neumann boundary conditions) corresponds to 9 = 0. Notice that there are only two 
boundary conditions, namely 9 = and 9 = oo, which are scale invariant. 

Let us now determine the operators T^ 1 (A) and K(\) defined in f|6. 2 j> . Since the de- 
ficiency indices of the operator A are n± = 1, the deficiency subspace Ker(^4^ — A) is 
generated by a normalized function which we denote by (j)\(x). Consequently, 

T^(\) : C -> Ker(^t - A) , 

(6.14) 

rr 1 (A)-c = A (x)/0 A (O)-c. 
The operator K(\) is therefore given by 

(6,5) ^ = -W>- 

The Krein's formula (Eq. ()6.7|) ) can then be written as 

{A - X)' 1 - (A°° - X)- 1 



.16) [A °-\) — (A°° — A) 



l + 0fT(A) 

Equation (|6.16|) gives the resolvent corresponding to an arbitrary selfadjoint extension 
A e in terms of the resolvents corresponding to the boundary conditions which are scale 
invariant, namely 9 = (Neumann) and 9 = oo (Dirichlet). 

Following [39], in the next section we will prove a generalization of the Krein's 
formula that relates in a similar way the resolvents corresponding to different self- 
adjoint extensions of the operator in (16.11) . The method employed follows the lines 
given by E. Mooers in |34j. We will obtain an expression similar to (I6.16P in which 
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the K(X) factor, although not given by (16. 15ft as in the regular case, is also related 
to the behavior near the origin of the functions in Ker (A^ — A) . 

6.2. Locally homogeneous second order differential operators. So, we con- 
sider the symmetric differential operator A 



(6.17) 

defined on C, 



A 



< + + V(x) 



C L 5 



). We assume that V(x) is an analytic function of 



x G R + bounded from below and the parameter v G (0, 1) C R. 

The following theorem describes the behavior of the functions in T>(A<) near the 
singular point x = 0: 

Theorem 6.4. Ifipe V(A^) then 

(6.18) ij(x) = C[$\ (x~ u+1/2 + 6$ x u+1/2 ) + 0(x 3/2 ) , 
for x — > + and some constants C[if)],9$ G C. 

Proof: By virtue of Riesz representation lemma 

(6.19) ip G V(A^) ->■ 3 i) G L 2 (R+) : (</>, A<j>) = $, 0) V0 G V(A) . 
Consequently, 

(6.20) A ] ip:=^p. 
If we define x := x~ u ~ 1 / 2 ip we obtain 

(6.21) d x {x 2u+1 d xX ) = -x u+1/2 (i> - V(x)iP) G Li(R + ) . 
Therefore, there exists a constant C\ G R such that 



(6.22) d xX = C x x- x - 2v - x- x - 2v 

The Cauchy-Schwartz inequality implies 



.1/+1/2 



, , v 2 - 1/4 



-9y + 



ipdy. 



-l-2u 



u+1/2 



(6.23) 

for some C 2 G 
(6.24) 



y 

<c 2 

In consequence, 



-9y + 



v 2 - 



111) 



tpdy 



< 



z/ 2 -l/4 



x 



(0,x) 



-l-2j/ 



.I/+1/2 



"3 + 



7/ 



l/4\ 



■0 ofa/ rfz 



<C 3 + C 4 x 



< 

1-1/ 
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where C3, C4 G M. Thus, there exist C^,Cq G R, such that 

(6.25) if) = C 5 a;"^ 1 / 2 + C 6 a^ +1/2 + 0(x 3/2 ) , 
for x — > + . 

□ 

Corollary 6.5. 

(6.26) M G P(At) (0, AV) - (AV,^) = C*[0]C[V] (9; - B*) . 



Remark 1: Notice that (16.261) verifies the first statement of Theorem 16.11 according 
to the definitions: Ti if) := C[if)] 9^ and T% if) '■= 1. 

Remark 2: By writing expression (16.261) for if) = we conclude that for all if) G 
P(A') the parameter ft^ defined by Theorem 16.21 is real. 

Proof: Expression 1 16. 26p follows from an integration by parts in its L.H.S. using 
expression 116. 181) . 

□ 

As a consequence of Corollary 16. 5[ the differential operator A admits a family of 
self adjoint extensions A e , characterized by the real parameter 9, whose domains are 
given by 

(6.27) V(A e ) := {0 G V(A^) : 9^ = 9} , 

where 9$ is defined according to Theorem 16.21 The parameter 9, with dimensions 
[length] ~ 2u , thus determines the boundary condition at the singularity. 

There exists another selfadjoint extension, which we denote by A°°, whose domain 
is given by, 

(6.28) D(A°°) = {(f) G V(A r ) : <j>(x) = C[<p] x u+1/2 + 0(x 3/2 ) , with C[<j>\ G C} . 

Let us point out that in the regular case limit, when v — > 1/2, where the sin- 
gular coefficient in the operator vanishes, this parameter 9 coincides with the one 
characterizing Robin boundary conditions for the regular case (see Eq. (I6.13P ). 

6.3. Generalization of the Krein's formula. Our purpose now is to establish a 
relation between the resolvents corresponding to different selfadjoint extensions of 
A. This relation will prove to be useful to show that the trace of the heat-kernel 
Tre~* A corresponding to a general selfadjoint extension admits, for 9 ^ 0, 00, a 
small-t asymptotic expansion with ^-dependent powers of t. 
We begin by stating the following theorem: 
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Theorem 6.6. For any f(x) G L 2 (M + ) and A ^ o~(A e ) there exists a unique function 
(f) e (x, A) G V(A e ) such that 

(6.29) (A e -\) ( j ) 9 (x,\) = f(x). 
Moreover, 

POO 

(6.30) <f> e (x,\)= I G e (x,x',X)f(x')dx' , 

Jo 

being Gg(x, x', A) the kernel of the resolvent {A 9 — A) -1 . 
The kernel Gg(x, x', A) can be written as 

Q(x' - x)L e {x, X)R(x', A) + Q(x - x')Lg(x', X)R(x, A) 



(6.31) G e (x,x',X) 



W[L e ,R](X) 



where B(-) is the Heaviside function. The functions Lg(x, A), R(x, A) satisfy equation 
(I6.29p for f(x) = 0. The latter is square integrable at x — > oo and the former satisfies 
the boundary condition 



3/2 ^ 



(6.32) L e (x, A) = x~ v+1 l 2 + 6 x v+l ' 2 + 0(x 

at x — > + . W[Lg, R](X) is the Wronskian of Lg(x, A) and R(x, A), and is independent 
of x. 

To obtain the generalization of the Krein's formula we begin by relating the re- 
solvents corresponding to 9 = oo and 8 = 0. In particular, for these selfadjoint 
extensions the boundary condition (16.321) reads 

(6.33) L^x, A) = x v+l/2 + 0(x 3/2 ) , 
and 

(6.34) L Q (x, A) = x- u+1/2 + 0(x 3/2 ) . 

Since these functions determine the behavior at the origin of the kernels Goo(x, x', A) 
and Gq(x, x', A), the following definitions are in order: 

Definition 6.7. 

(6.35) GocO^A) ■=limx- v - 1/2 G 00 (x,x',X) , 

x—>-0 

(6.36) G (x', A) := lim x u ' 1/2 G (x, x' , A) . 

These new functions Goo(x, A) and Go(x, A) determine the behavior at the singu- 
larity of the solutions 0°°(x, A) and 0°(x, A) of (16.291) corresponding to 9 = oo and 
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= 0, respectively. Indeed, 

POO 

(6.37) 0°°(a;,A) = / Goo(x,x', X)f(x') dx' = 0°°(A) x v+1/2 + 0(x 3/2 ) , 



(6.38) 4>°(x,\) 

being 

(6.39) 

(6.40) 



G (x, x', X)f(x') dx' = (f>°{\) x~ v+1/2 + 0{x 3 / 2 ) 



<F(A) := / G OQ (x',X)f(x')dx' 



0°(A) := / G (x',X)f(x')dx' . 
'o 



To obtain a relationship between the kernels G^x, x', A) and G (x, x', A) we will 
relate the solutions (p°°(x, A) and <f)°(x, A) of (16.291) corresponding to the same inho- 
mogeneity f{x). To do that we need the following two lemmas: 

Lemma 6.8. Let ip (x) E V(A°) such that <p (x) = x~ v+l ' 2 + 0{x 3/2 ) for x -> + . 

Then, the solutions 0°°(x, A) and 4>°(x, A) of K6. 29\) are related by 

(6.41) 



<r(x,A) = 0°(x,A)-0°(A) 



<Po(x)- I G OQ (x, x', X)(A° - X)(p (x') dx' 
'o 



Proof: On the one hand, 

(6.42) (A* — A) (p°°(x, A) — f(x). 

Moreover, 

(A^ - A) |</>°(o;, A) - 0°(A) (p (x) - G^x^' , X)(A° - X) Vo (x') dx' 

= f(x) - 0°(A) [(A - AV (x) - (A - \)<p (z)] = f(x) . 

(6.43) 

On the other hand, for x — > + , 

(6.44) 0°°(x, A) = 0°°(A) x"+ 1/2 + 0(x 3/2 ) , 

and 



0°(x,A)-0°(A) 
= 0°(A) 



<^o(^) - / Goo(x,x', X)(A° - X)ifo(x') dx' 



GnixWKA -X)<p (x')dx' 



x u+l/2 + Q ^ x 



3/2 \ 



(6.45) 
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As a consequence, both sides of equation (I6.4ip belong to V(A°°) and satisfy Eq. 
(I6.29P for 6 = oo. Expression (16.411) is then proved by virtue of the uniqueness 
stated in Theorem 16.61 

□ 

Lemma 6.9. 

POO 

(6.46) <p (x)- / G x (x,a/,\)(A -\)(p (a/)<y = 2vG oo (x,\). 

Jo 

Proof: Since the kernels of the resolvents G 00 (x, x', A), Gq(x, x', A) are symmetric 
(see Eq. (IQTl) ) 

(6.47) A* [G (x, x', A) - G^x, x', A)] = , 
with A* acting either on x or x' . Thus, 

poo 

<Po(x) — / G oo (x,x',X)(A - X)(p (x')dx' = 



[G (x, x', A) - Goo(x, x', A)] (A - A)^o(s') <te' = 
= lim {[G (x,x',X) — G 00 (x,x',X)\ ■ ip (x)- 

x— H)+ 

- d x > [G (x, x', A) - Goo(x, x', A)] • v?oO)} = 
= [x~ v+l/2 G (x, A) - x v+l ' 2 GooCx, A)] ■ (-1/ + 1/2) aT^ 1 / 2 - 
- [(-x/ + 1/2) x-^ 1/2 G (x, A) - (u + 1/2) x u ~ 1/2 G^x, A)] • aT^ 1 / 2 = 
(6.48) = 2vG 00 (x,X). 



□ 



Lemmas 16.81 and IB31 lead to the following result: 
Lemma 6.10. 

(6.49) <P°{x, A) = <p°°{x, A) + 2v G^x, A)0°(A) . 

We are interested in rewriting expression ACT]) so that (j)°(x, A) be given in terms 
of quantities corresponding to the extension characterized by 9 = oo. In so doing, 
we take the x — > + limit in equation (I6.49P obtaining 

(6.50) G^x, A) = — (x- u+1/2 - imrV +1/2 ) + 0(x i/2 ) , 

2v 

where 

(6.51) K(X) r- /(A) 



0°°(A) ' 

The term K(X) defined in (I6.5ip relates the behavior at the singularity of the 
solutions to equation f !6.29|) corresponding to the selfadjoint extensions 9 = oo and 
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= 0. Notice that Eq. (I6.50p allows us to compute K(X) by studying the behavior at 
the singularity of the kernel of the resolvent corresponding to the extension 9 = oo. 
Therefore, the kernel Goo(x,x' , X) determines K(X) and, consequently, also <f> (A). 
Notice that K(X) reduces to expression (16.151) for the regular limit v — > 1/2. 

We can finally express the solution A) to (I6.29j) corresponding to 9 = by 
means of the data obtained by imposing the boundary conditions corresponding to 
9 = oo (see Lema |6~T0]) . 

(6.52) <f>°(x, X) = 0°°(x, A) + 2uK(X) C7 00 (x, A)0°°(A) . 

Since this equation is valid for any inhomogeneity f(x), by virtue of equations (16.371) . 
(I6.38P and (I6.39p . we obtain the following theorem: 

Theorem 6.11. 

(6.53) Go(x, x , A) = Goo(^, x', A) + 2uK(X) Goo(x, A)C7 00 (x / , A) . 

Next, we will establish an expression similar to (I6.53j) giving the resolvent for an 
arbitrary selfadjoint extension in terms of data related to the boundary conditions 
corresponding to 9 = oo. The first step is to state the following lemma: 

Lemma 6.12. The solution <f) e (x,X) to Ii6. 29\) is given by 

(6.54) <f) e {x, A) = 0°°(x, A) + 2u (K(X)- 1 + 9)' 1 C7 00 (x, A)0°°(A) . 

Proof: By means of equation (I6.49P it is immediate to show that the difference 
between both sides of expression (I6.54p belongs to Ker(y4^ — A). 

On the other hand, both sides of (I6.54p belong to T>(A d ) since the behavior of its 
R.H.S. at the singularity is given by (see eqs. (16.371) . (I6.50P and (16.511) ) 

Once more, uniqueness established in Theorem 16.61 leads us to equation (16.541) . 

□ 

From Lemma I6.12[ together with equations (16.301) , (I6.37P and (I6.39P , we straight- 
forwardly get the following theorem: 

Theorem 6.13 (Generalization of Krein's formula's). 

(6.56) G e (x } x', A) = G^x, x', A) + 2u (K(X)- 1 + 9)' 1 C7 oc (x, X)G 00 (x', A) . 
Expressions (I6.53P and (I6.56P readily lead to 

1 r nn\ n / I \\ n I I \\ % > A) G^ (x, X , A) 

(6.57) G e (x,x , A) - Goo(x,x , A) = \ + 9K(\) ' 
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Therefore we obtain the following relation between the resolvents of the different 
selfadjoint extensions: 

( , 58) v-ir-^- v-T^-** - 

This expression formally coincides with the Krein's formula (I6.16p . which is valid 
for regular operators. However, while the factor K(X) in (I6.16P is given by (16.151) . 
in the singular case under study K(X) in (I6.58P corresponds to equation (I6.5ip . As 
already mentioned (I6.16P and (I6.58P coincide in v — > 1/2 limit. 

We resume our results in the following theorem, which will allow us to prove the 
non-standard behavior of the spectral functions of the SAE of the operator A in 
fl6TTD: 

Theorem 6.14. 

Tr{(A°- A)- 1 - (A°° - Xy 1 } 



(6.59) Tr {(A 9 - X)- 1 - (A°° - X)- 1 } 



1 + 6K(X) 



In the next section we will show that the asymptotic expansion of K(X) for large 
|A| presents powers of A whose exponents depend on the parameter v. This will 
finally lead to the asymptotic series (16.21) . 

6.4. Asymptotic expansion of the resolvent. In this section we will make use 
of Theorem 16. 141 to obtain the large- |A| asymptotic expansion for the resolvent (A — 
A) -1 of an arbitrary selfadjoint extensions of the operator A defined in Eq. (16.171) . 
According to this theorem it suffices to study the solutions to 

(6.60) (A + z)ip = 0, 

satisfying the boundary conditions corresponding to 6 = oo and 6 = 0. If we consider 
A in the negative real semi-axis we can take z G M + . In particular, we will focuss 
on the behavior of the solutions for large z. 

Taking into account the scaling properties of the first two terms in (16.171) it will 
be convenient to define a new variable y := yfzx G M + . The solution to equation 
(I6.60p can then be written as if) = ip(y/zx, z), being ip(y/zx, z) a solution to 

(6.61) (-dl + + 1 + \ v{y/V~z)\ Hv, z) = o. 

We propose the following Ansatz 

oo 

(6.62) if;(y, z) = <p(y) + ^ ^ n {y)z~ x - n l 2 , 

n=0 
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(6.63) 



x 



n=0 



where V n := V^(Q)/n\. Eq. flTOD can be now solved order by order in z. The 
solution to (I6.6ip which is square integrable at y — > oo can be written as 



(6.64) 

where 
(6.65) 



R(y, z) = <JyK v {y) + J^Mv)^ 



/2 



n=0 



- y 1/2 h{y) f 



4>n(y) = 

ftflvirf) dy' - 
y'K v {y') dy' 



V n y'^^K v {y') + Viv'^M) 

l+m=n—2 

V n y^ +1 ^K v {y') + Vw'^mW: 

l+m=n—2 

Therefore, the behavior of R(y, z) at y — > + is given by 

(6.66) R(y, z)c^± y^ + 2=£ H(z) • y^ + ... 

where we have defined 
(6.67) 



TT , . , 2sin(7rz/) 
H(z) := H — - x 



7T 



oo p C 



V n y n+1 ' 2 K u {y) + Vw^Uv) 



l+m=n—2 



VyK u (y) dy. 



It is important to notice that H(z) admits a large- z asymptotic expansion in half- 
integer powers of z. 

Next, we find a relation between K(z) in equation (I6.59j) (defined in Eq. (16.5ip ) 
and H(z). To obtain an expression for K(z) we study the kernel of the resolvent 
Goo(x, x',z) which, for x < x', is given by (see Eq. (I6.3ip ). 



(6.68) 



-1/2 



W[L,R}(z] 



-L(y,z)R(y',z) 



y=y/~zx,y' =y/zx' 



The function L(y, z) is a solution to ( I6.6ip whose leading term at the origin is 
proportional to y u+1 ^ 2 . W[L, R](z) is the Wronskian of L(y,z) and R(y,z), which 
is independent of y. 
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According to definition ( 16. 35ft 
(6.69) 



OQ \*£ 1 



W[L,R](z) 



L(y,z)R(y',z) 



y=0,y'=^/zx' 



Replacing (I6.66j) into equation (16.691) we obtain the behavior of G^x, z) for x — > + , 



(6.70) 



G^(x ~ 0,z) 



z -^ y ---^L(y,z) 



y=0 



^l(Vzx)-^ + H{z) ■ (V~z^ +1/2 



W[L,R](z) 
^zx) 



+ . . 



2 1 ~ v ^ ' 2 1 
Comparing equations (16.501) and (I6.70p we get a relation between K(z) and H(z) 



(6.71) 



Since H(z) admits an asymptotic expansion in half-integer powers of z, the large- z 
asymptotic expansion of K(z) contains powers of z whose exponents depend on the 
parameter v. Theorem 16.141 shows that this powers are also present in the large- 
z asymptotic expansion of the resolvent trace, which after equation (I6.59P can be 
written as 

Ti{{A° + z)- 1 - {A°° + z)- 1 } 



(6.72) 



Tr{(A e + z)- 1 - (A^ + z)- 1 



} 



The trace in the R.H.S. of equation (I6.72p can be readily obtained from equation 
0633]) 

poo 

(6.73) Tr{(A° + z)' 1 - (A°° + z)- 1 } = 2uK(z) G 2 oc {x,z)dx. 

Jo 

To evaluate this expression we compare equations (I6.50P and ( I6.70P and we get 

z -l/2 y -v-l/2 L fa ^ 

< 674) WKWT 

Therefore, equation ( 16.691) reads 



1 



2 v uT(v] 



u-l/2 



(6.75) 



Goo yl' 5 Z) 



■u-l/2 



2 lJ vY{v) 

Replacing this equation, together with ( 16. 71 p . into expression (16. 73 ft we obtain 
(6.76) Tr {(A + z)-' - (A™ + z)-'} = — U- -/ R^zfdx. 

r>)r(i - v) j 
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This shows that Tr {(A + z)^ 1 — (A 00 + z)" 1 } admits an asymptotic expansion in 
half- integer powers of z. 

The following theorem summarizes equations (16.721) and (I6.76P regarding the 
resolvent-trace. We also state the corresponding result as regards the heat-kernel 
trace, which is the inverse Laplace transform of the resolvent-trace. 

Theorem 6.15. 



• The trace of the difference between the resolvents (A e — A) 1 and (A 00 — A) 1 
admits an asymptotic expansion for large |A| given by, 

oo 

Tr{(A e - A)" 1 - - A)" 1 } ~ V ) A ~ f + 

n=2 

oo 

(6.77) + 0N,n(v, V) 9 N A -1 



2 2 

N,n=l 



The coefficients a n (i>, V), (3 n (i>, V) depend on the parameter v characterizing 
the singularity and are also determined by the coefficients V n characterizing 
the analytic potential V(x) by means of equations \6. 72$ , Ij6. 76$ , &6.64$ , 
(ESf and (fOTP with z = e™\. 



• The trace of the difference e tA — e tA °° admits a small-t asymptotic expan- 
sion given by, 

(6.78) 

1 ] Z^l>iV+f + §) 

□ 

Let us give a dimensional analysis argument to explain the non-standard powers 
of t in the expansion (16.78 p . First of all, notice that the parameter 9 introduced 
by the boundary conditions has dimensions [length] ~ 2u and, after the analyticity of 
V(x), the dimensions of every other parameter in the problem is an integer power 
of the length. Since t has dimensions [length] 2 , if the coefficients of the asymptotic 
expansion of the heat-trace were to depend analytically on 9, then it is necessary 
that this expansion contains integer powers of t u . The only selfadjoint extensions 
for which these powers are to be absent are those with 9 = and 9 = oo. 

Example: Let us consider V{x) = x 2 . We will use the expansion (I6.78P to describe 
the pole structure of the difference between the ^-functions Ca( s ) — C^°( s ) corre- 
sponding to the operator (I6.17p . which will confirm the results obtained in Section 
[5] with two other techniques. 
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First of all, notice that for V(x) = x 2 only one of the coefficients V n defined in 
( I6.63P is non vanishing, namely 

(6.79) V n = 5 n , 2 . 

As a consequence, the only functions if> n {y) which are non-trivial corresponds to 
n = 2 + 4k with k = 0, 1, . . . (see Eq. (I6.65P ). According to equations ( I6.67P and 



(6.80) 



(6.8i; 



oo 

k=l 

oo 



2k 



-2k 



k=l 



for some Ck(is),C k (u,y). Substituting these equations into ( 16.761) we get 



(6.82) Tr{(A° 



(^°° + z)- 1 } ~ v z- 1 + J2 C '^ 



V) z 



-3-2fc 



k=0 



where C k (u) can be written in terms of Ck{y),C' k {y,y). 

Replacing equations ( 16. 82ft and ( 16.801) into ( 16.721) we obtain the asymptotic expan- 
sion of the trace of the difference between the resolvents {A 6 + z) -1 and (A°° + z)~ l 



(6.83) Tr { (A e + z)' 1 - (A°° + z)' 1 } 



v z 



V) z 



-3-2k 



k=0 



X 



N=0 



ra-v) 



-i N 



N -Nu 



e N z 



kW) Z 



-2k 



k=l 



N 



It is straightforward to see that the first series in ( 16.83}) gives no contribution to the 
pole structure of the difference Ca( s ) ~ (a( s )- O n the other hand, from the second 
series in the asymptotic expansion (I6.83j) one shows that Ca( s ) ~ Ca( s ) h &s simple 
poles which are located at 

(6.84) s Nin = -Nu - 2n with N — 1,2, . . . and n = 0, 1, . . . 

This result coincides with Eq. (I5.27[) . In particular, the leading term in (I6.83|) leads 
to a simple pole at 

(6.85) s lfi = -u, 
whose residue is given by 

(6-86) f^O, 
in agreement with the result quoted in Eq. ( I5.28[) . 
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7. Conclusions 

We have studied some symmetric non essentially self adjoint first and second order 
differential operators with a singular potential term with the same scaling dimension 
as the highest derivative term, characteristic which we have mentioned as local scale 
invariance (at the singular point in the potential). For certain range of the "external 
parameters" that weights the singular term (i.e., the coupling constants in the po- 
tential - g or v throughout the paper) these operators admit a continuous family of 
selfadjoint extensions. Each selfadjoint extension describes a different physical sys- 
tem, with its own spectrum determined by the behavior of the functions belonging 
to its domain near the singularity. This is an essential point in our discussion for it 
implies the existence of infinitely many physically admissible boundary conditions at 
the singularity, identified in our examples by an additional dimensionful parameter 
(which is not present in the expression of the differential operator). 

This situation takes place when the adjoint of the differential operators considered 
above admits two different square integrable behaviors near the singularity for the 
functions in its domain of definition. Then, a symmetric extension will contain, in 
general, both possible behaviors and then it must incorporate an additional dimen- 
sionful parameter necessary to specify its domain. The only exceptions are those 
combinations which make the domain (locally) scale invariant. 

Therefore, the spectral functions associated to a selfadjoint extension of these 
operators will, in general, depend on this additional dimensionful parameter (with a 
scaling dimension which depends on the external parameters in the singular potential 
term), opening the possibility of having, for example, non-standard powers of t in the 
asymptotic expansion of the the heat-kernel trace and, consequently, non-standard 
poles in the associated ^-function. By non-standard we mean that these powers and 
poles are not determined by the dimension of the base manifold and the order of 
the differential operator only (as in the smooth coefficients case) but also depend 
on the external parameters that characterize the singularity. 

Our results show that those selfadjoint extensions which break in this way the 
local scale invariance of the dominant scaling dimension terms in the differential 
operator present these non-standard poles in the associated (^-function as well as 
non-standard powers of A in the large- 1 A | asymptotic expansion of the resolvent. As 
a consequence, for second order differential operators one finds non-standard powers 
of t in the small t-asymptotic expansion of the heat-kernel trace. Moreover, the 
Seeley-De Witt coefficients and the residues of the (- and rj- functions depend on 
the selfadjoint extension. 

In fact, also the large- n asymptotic behavior of the eigenvalues X n of the selfad- 
joint extensions contains powers of n which depend on these external parameters. 
One also expects non-standard singularities in the corresponding Green functions at 
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coincident points. This issue is relevant for the definition of physical states and the 
regularization of the stress-tensor, for example. 

These non-standard behavior of the spectral functions has been explicitly shown 
by solving some examples of first and second order non essentially selfadjoint differ- 
ential operators on both compact and non-compact one-dimensional base manifold. 
We have also proved that this phenomenon is not affected by the introduction of 
arbitrary smooth potentials. 

Let us finally mention that, in establishing these results, we have derived an exten- 
sion of the Krein's formula which applies to this kind of differential operators defined 
on functions which, generically, do not have a regular behavior at the singularity. 
This formula (Eqs. ( I6.56P and ( I6.58P ) relates the resolvents of different selfadjoint 
extensions and directly leads to the non-standard small-t asymptotic expansion of 
the heat-kernel trace we have found. 



Acknowledgements: We are very pleased and we feel very honored for the pos- 
sibility of contributing to this special volume in honor of Prof. Stuart Dowker. 
This work was supported in part by grants from CONICET (PIP 01787), ANPCyT 
(PICT 00909) and UNLP (Proy. 11/X492), Argentina. 



Given an elliptic differential operator A in a manifold M with a complete or- 
thogonal set of eigenf unctions corresponding to the eigenvalues {A n } nS N, then the 
associated ^-function is defined as [9j [TO] 



rteN 

which is a convergent series for TZ(s) large enough. 

If |A n | — > oo fast enough when n — > oo, then also exists the trace of the resolvent, 
given by 



Moreover, if the set {9ft (A n )} is bounded below and 9ft (A n ) — > oo for n — >• oo, then 
the trace of the associated heat-kernel is given by 



Appendix A. Spectral functions and its relations 



(A.l) 



( a (s)=TtA- s = J2K s , 



(A.2) 




(A.3) 




neN 



For positive definite operators we have for the Laplace transform 



(A.4) 
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for TZ(X) < X n ,Vn, and for the Mellin transform 



OG 



(A.5) TrA~ s = -^-l t s ~ l Tre- tA dt, 

r(s) Jo 

for TZ(s) > m/d, where m is the dimension of the manifold and d the order of the 
differential operator. 

The (^-function singularities are related with the asymptotic expansion of Tr (A — 
A) -1 for large |A| and with the asymptotic expansion of Tre~ tA for small values of 
t. Indeed, if Ca( s ) has simple poles at s = s n < so, for n G N, then Tr (A — A) -1 
has an asymptotic expansion in powers of the form A Sn_1 , while Tre~ tA admits an 
asymptotic expansion in powers of the form t~ Sn . In particular, under the hypothesis 
considered in [9], these powers depend only on d and m. Moreover, the coefficients 
of both expansions are determined by the residues of Ca( s ) & t the corresponding 
poles. 

For example, if 



(A.6) Tr{e- tA }~J2 c ^ f 



n=0 

with s n < So, then for 3f?(s) > Sq we have 

n 



r( S )( A (s) = J2cn(A) / + 
n=0 J( > 

t 3 - 1 f Ti{e- tA } - Y^c n {A)r s A dt+ 



no 

h(s) 



k fo \ S - S nJ 

where h(s) is analytic on the open half plane 3ft(s) > s no . Therefore, the residue of 
r(s)Ci( s ) a t s n is related with the coefficient c n (A) through the equality 

(A.8) Res[T(s)( A (s)}\ s=Sn = c n (A). 

Since T(s) has simple poles at s = 0, —1, —2, . . . , the residues Res [Ca(s)]| s vanish 
when s n = 0, — 1, —2, .... In particular, (a{s) is analytic in a neighborhood of the 
origin. 

Even for nonpositive elliptic differential operators A, the complex s-power of A is 
defined in terms of the resolvent as [9] 

(A.9) A' 3 = — l\- s {A- Xy 1 dX, 

2 m Jr 
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where C is a curve enclosing anti-clockwise the eigenvalues of A. From this one gets 

(A.10) TiA- s = — (f\- s Ti(A- Xy 1 dX. 

2 it i J c 

The resolvent (A — A) -1 , the complex power A~ s and the heat-kernel e~ tA can be 
considered as integral operators characterized by its kernels, G(x, x', A), (a( x , x ' , s ) 
and K(x, x', t), defined for A G C\{A n } n£ N, 72.(s) sufficiently large and t > respec- 
tively. In this case, their traces are expressed as 

(A.ll) Tr(A-X)- 1 = f G(x, x, A) dx , 

Jm 

(A.12) Ca(s) := Tr A~ s = / ( A (x,x,s)dx 

Jm 

and 

(A. 13) Tre~' A = / K(x, x, t) dx . 

Jm 



Appendix B. The Hankel expansion 

In this appendix we write down some of the asymptotic expansions of the Bessel 
functions that are used in Section ETT1 [55]. For \z\ — > oo, with v fixed and | argz| < ir, 
we have 

(B.l) J v (z) ~ ( — ] {P(v, z) cos x(v, z) - Q(v,z) sin x(v,z)} , 



IT Z 



where 



(B.2) x{v,z) = z-fa + \)*, 

oo 

(B.3) P(z,,^)~^(-l) fc (z,,2A;) — 



fc=0 

oo 



(2z) 2k ' 



(B.4) Q(„, z) ~ £(-1)*<i/, 2k + 1) —L^ 



and the coefficients 



(B.5) (u,k)= ^^^ + k \ =(-u,k) 

ft! T (J + v - k) 
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are the Hankel symbols. Therefore, P(—v,z) = P(v,z), Q(—u,z) = Q(v,z) and 



(B.6) (|- 

where the upper (lower) signs correspond to z in the upper (lower) open half-plane. 
In particular, the quotient 

(B.7) ~ g ; ' ~ e ±m U-s) , 

for Of (A) > and Q'(A) < 0, respectively. 

Similarly, the derivative of the Bessel function has the following asymptotic ex- 
pansion for | argz| < n, 

(B.8) Jl(z) ~ --^={R(v, z ) sin X (v,z) + S (is, z) cos X (v,z)} , 

where 

^ x , v z/ 2 + (2A;) 2 - 1/4 (i/,2fc) 

(R9) ^^g(-l) ^ 2 _^_ v / )2 ^, 



(B - 10) ^' " } , 2 -(2A; + l-l/2) 2 • 

Then, 



k=0 



(B.ll) J' v {z) ~ =r = z) =f i 5(1/, z)} , 

V27T2; 

where the upper sign is valid for $f(z) > and the lower sign for $$(z) < 0. From 
the following relation 

(B.12) R(v,z)±iS{v,z) = P(is,z)±iQ(is,z)+T±(is,z) , 

with 

k 



(B.13) T ± M~£(2A;-1)<^-1) ^ , 

fc=i ^ ' 



we get 



2 
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where the upper sign is valid for ^s(z) > and the lower one for ^s(z) < 0. Finally, 
since the Hankel symbols are even in v (see Eq. flB.5j) ) we have 

(B.i5) m„*M. 

Appendix C. Closure of H 

In Section [5] we omit the contributions to the boundary condition in Eq. (I5.12p of 
the functions in the domain of the closure of the operator H defined in Eq. (j5.ip ; in 
this appendix we will justify this procedure. Indeed, we will show that if e D(H) 
then 

(C.l) <j>{x) = o{x a ) and <j>'{x) = o{x a ^) 

near the origin, for any a:=z/+i<3/2. 

In order to determine the domain T>(H) of the closure of H we must consider 
those Cauchy sequences {(p^nen in T^{H) = C^°(1R + ), such that {H(p n } n£ ^ are also 
Cauchy sequences. Since the coefficients of H are real, it will suffice to consider real 
functions. Throughout this section ip := ip n — (p m (with n, m G N) so that ip — > 
and Hip — > as n, m — » oo. 

In the following we write g := z/ 2 — ~. Notice first that the scalar product 

(C.2) (<p,H<p) = J (V 2 + J^ 2 + xV) dx< \\<p\\\\H<p\\->0 

as n,m — > oo; || • || represents the usual norm in L2(M + ). Therefore, for g > 0, we 
conclude that 

(C.3) {ip' n (x)} neN , y n X i and (x)} neN 

are also Cauchy sequences. We will now prove the following 

Lemma C.l. Let {<^„} n eN be a Cauchy sequence in T>(H) = Qj°(R + ) such that, for 
g > 0, 1 < a < 2 and g ^ (a 2 - l)/4, 

(C.4) {^„}„ 6N , (^Mj , and W W 



r.a— 1 



are a/so Cauchy sequences. Then, 
(C.5) I^U and^»W 



rpl+a/2 \ 1 rv.a/2 



are Cauchy sequences, too. 
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Proof: Taking into account that the sum of fundamental sequences is also a Cauchy 
sequence, we see that 



(C.6) 



as n,m — >■ oo, for any pair of real numbers A and B. It is easily seen that appro- 
priately cho 
the lemma. 



priately choosing the coefficients A,B (whenever g = (a 2 — l)/4) Eq. ( 1C6|) proves 



□ 

Let us now assume that g is an irrational number. Then, applying iteratively 
Lemma fC. II to the sequences (10.3|) one can show that, for any positive integer k, 



(C.7) 



2[l-(l/2) fc ] 



X 



and 



neN 



„2[l-(l/2) fc ]-l 



neN 



are Cauchy sequences, too. One immediately concludes that {x +E cp n (x)} nen and 
{x~ l+e <p' n (x)} nen are also Cauchy sequences. If g were a rational number, one could 
choose an irrational a G (1,3/2) from which Lemma [C.ll could also be iteratively 
applied to arrive to the same conclusions. 

In the following we will consider the behavior of the functions near the origin. For 
any e > 0, we can write 



(C.8) x->(x)=/ {y- a v(y))'dy= y~ a+1 



-a 



<p(y) , <p'(y) 



,2-e 



+ 



,1-e 



dy. 



Therefore, for x < 1, a < 3/2 and e small enough, we have 



(C.9) 



\x a (p(x) I < 



1 \ 1/2 

y 2(-a+l~e) dy 







a 



<p(y) 



y 



2-s 



+ 



y 



l-e 







as n,m — > oo. We conclude that the sequence {x~ a v? n (x)}„ e N, for a < 3/2, is 
uniformly convergent in [0, 1] and its limit is a continuous function vanishing at the 
origin, which we write as x~ a 4>(x): 



(CIO) 
(C.ll) 



lim (x a ip n (x)) 



x 



[X 



lim (x~ a 6(x)) = . 



In particular, for a = we have the uniform limit 
(C.12) lim <p n (x) = (j>{x) , 

n— >oo 

which coincides with the limit of this sequence in L 2 ( 
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Similarly, we can write 
(C.13) 



x- a+1 ip'(x) = - 



-a+l 



Hip(y) dy + 



+ 



y 



-a+l-e 







- a+l) ^ +g ^\ dy+ h 



y a+2 yp(y)dy. 



Therefore, for x < 1, a < 3/2 and e sufficiently small, we have 

•1 \ 1/2 



(C.14) 



|aT a+ V(aO| < 

1 \ 1/2 

y 2(-a+l-e) dy 



y 



la — 11 



2(-a+l) 



dy) ||^(y)|| + 



y 



2(-a+2) 



,1-e 



+ 3 



x 



2-e 



+ 



1/2 



as n, m— t-oo. Consequently, the sequence ^(x)}^^, with a < 3/2, is 

uniformly convergent in [0, 1] and its limit is a continuous function vanishing at the 
origin, which we write as x~ a+1 x( x ) : 



(C.15) 
(C.16) 



lim (x- Q+1 <p' n (x)) =x- a+1 X (x), 

n— »oo 

lim (aT Q+1 x(x)) = . 



In particular, for a = 1 we have the uniform limit 
(C.17) lim <p' n {x) = X {x), 

n— ¥oo 

which coincides with the limit of this sequence in L2(1R + ). Let us now show that 
x(x) = 4>'(x). Indeed, for x < 1, we have 



(C.l* 



< \4>{x) - (f n (x) \ + 



[X - 



< 



(x(y) - <p'M) d y 



x(y) dy 

< \(f)(x) - <p n (x)\ + \\x - (p' n \\ -»> o 



as n — > oo. Then, 4>(x) is a differentiable function whose first derivative is x( x )- On 
the other hand, eqs. ( 1C.11I) and ( 1C.16I) imply that any e T^{H) satisfies Eq. ( 1C.1I) . 
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